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Chapter  1 
INTRODUCTION 

1.1   Discussion  of  problem 

Water  is  one  of  the  materials  required  to  sustain  life 
and  has  long  been  suspected  of  being  the  source  of  many  of 
the  illnesses  of  man.  It  was  not  until  a  little  over  100 
years  ago  that  definite  proof  of  disease  transmission  through 
water  was  established.  For  many  years  following  the  maior 
consideration  was  to  produce  adequate  supplies  that  were 
hygenically  safe.  The  public  has  been  more  exacting  in  its 
demands  as  time  has  passed  and  today  water  engineers  are 
expected  to  produce  finished  waters  that  are  free  of  color, 
turbidity,  taste,  odor,  harmful  metal  ions  and  organic 
compounds. 

It  has  long  been  known  that  all  natural  bodies  of  water 
have  the  ability  to  oxidize  organic  matter,  provided  that  the 
organic  loading  rate  is  kept  within  the  limits  of  oxygen 
resources  in  water.  Owing  to  the  wide  use  of  pesticides  and 
the  use  of  chemicals  by  various  industries  there  is  a  great 
deal  of  organic  contamination  in  surface  and  ground  waters. 
Most  of  these  organic  compounds  have  adverse  influence  on 
aquatic  and  human  life. 

Every  organic  contaminant  cannot  be  regulated .  Hence  the 
EPA   /  Environmental  Protection  Agency  /  decided  to   regulate 
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2 

these  compounds  based  on  the  following  considerations : - 

1)  Physical  volume  produced  yearly 

2)  Physical  and  chemical  properties 

3)  Estimated  human  exposure 

4)  Toxicological  suspicion  and  opinion 

5)  Public  interest 

6)  Significance  to  society 

Based  on  the  above  criteria  EPA  issued  a  list  of  pollutants 
which  are  called  priority  pollutants.  In  this  list  there  are 
128  compounds  of  which  14  are  metals  and  114  are  organic 
compounds.  Of  these  114  organic  compounds  11  are  phenolic 
compounds. 

Phenolic  compounds  are  of  growing  concern  as  water 
pollutants.  Some  phenols  have  very  low  taste  and  odor 
thresholds  (eg  .2  ^cg/L  for  2-chlorophenol  )  whereas  others 
are  highly  persistent  and  toxic  eg .  pentachioro- phenol  <PCF) 
Buikema  et.  al  (1979)  reviewed  the  occurence  and  biological 
effects  of  phenolic  compounds  in  aquatic  ecosystems.  More 
than  90&  of  phenols  are  produced  by  industrial  synthesis „  but 
some  phenolic  compounds  are  natural  products.  Phenols  are 
used  for  example,  for  wood  preserving  or  as  chemical 
intermediates  for  dye  manufacturing.  Natural  phenols  are 
formed  in  aquatic  and  terrestrial  vegetation  and  can  be 
released  as  pollutants  by  the  pulp  and  paper  industry. 

Bueikma  et.al  C  1979  ) ,  Van  Hall  C  1979  )  have  studied 
the  occurences  of  phenolic  chemicals  in  various  tyoes  of 
industrial   wastewaters.    Hexachlorophene  and  PCP  have   been 


investigated  in  sewage  and  water  by  Butler  et  al  .  (1330). 
Dietz  and  Trand  have  determined  large  number  of  phenols, 
especially  chlorinated  derivatives  in  domestic  sewage 
effluents,  various  wastewaters  and  surface  waters. 

Chlorinated  phenols,  otherwise  known  as  chlorophenols, 
are  readily  soluble  in  water.  Most  of  these  chlorophenols  are 
carcinogenic  in  nature  and  they  produce  sicknesses  such  as 
depression  of  the  central  nervous  system,  increase  in  weight 
of  the  liver  and  kidney,  gastrointestinal  problems  etc.  Some 
of  the  industries  that  have  chlorophenols  in  their 
wastewaters  are  textile,  ink,  metal  finishing,  steam 
electric,  leather  tanning  and  finishing,  petroleum  and  paint 
industries. 

Owing  to  the  health  problems  exhibited  by  organic 
compounds  on  life  and  since  they  impart  taste  and  odor  to 
waters,  their  removal  from  waters  and  wastewaters  is  a  must. 

1 . 2   Methods  of  Removal 

There  are  two  methods  by  means  of  which  organic 
compounds  are  removed  from  waters  :  1>  Aeration,  and  2) 
Adsorption. 

Aeration  or  adsorption  can  be  used  to  remove  volatile  organic 
compounds,  but  organic  compounds  with  high  solubility  can  be 
removed  from  waters  by  adsorption  quite  effectively, 

1.2.1  Aeration 

Aeration   is  usually  done  by  passing  air  under   oressure 


through  a  body  of  water.  Aa  the  air  bubbles  up  to  the  surface 
it  removes  the  volatile  organic  compounds  present  in  the 
water  and  discharges  it  into  the  atmosphere.  However  with 
aeration  it  is  not  possible  to  remove  volatile  organic 
compounds  present  in  water  to  low  concentrations. 

1.2.2  Adsorption 

Adsorption  is  defined  as  the  tendency  exhibited  by  all 
solids  to  condense  upon  their  surface  a  layer  of  any  gas  or 
liquid  with  which  they  are  in  contact.  Activated  carbon  is 
used  extensively  for  adsorptive  purposes  because  of  its 
termendous  surface  area  in  relation  to  its  mass.  The  total 
surface  area  represented  by  the  walis  of  its  pores  is  in  the 
range  of  900-1300  square  meters  per  gram  of  carbon. 

Nearly  every  organic  material  capable  of  being 
carbonized  is  a  possible  candidate  for  the  production  of 
activated  carbon.  Among  the  most  common  are  wood.  coal, 
nutshells  and  petroleum  residues. 

Activation 

Activation  is  a  physical  change  wherein  the  surface  of 
the  carbon  is  tremendously  increased  by  the  removal  of 
hydrocarbons  from  the  carbon  structure.  There  are  several 
methods  available  for  the  activation  of  charcoal.  The  most 
widely  employed  are  treatment  of  the  carbonaceous  material 
with  oxidizing  and  purging  gases  such  as  air,  steam  or  carbon 
dioxide  and  the  carbonization  of  raw  material  in  the  Dresence 


of  chemical  agents  such  as  zinc  chloride  or  phosphoric  acid. 

1.3  Parameters  that  influence  Adsorption 

There  are  several  parameters  which  influence 
adsorption : - 

1)  Temperature 

2)  pH 

3)  Nature  of  adsorbate 

4)  Particle  size 

Generally,  adsorption  capacity  increases  when  temperature 
decreases  or  when  pH  increases.  Adsorption  depends  also  on 
the  molecular  weight,  solubility,  and  the  concentration  of 
the  adsorbate.  As  the  adsorbent  particle  size  decreases  the. 
rate  of  adsorption  is  generally  higher. 

Granular  activated  carbon  adsorption  is  an  effective 
unit  process  for  the  removal  of  synthetic  organic  chemicals, 
trihalomethanes  and  other  organic  haiides  from  drinking 
water.  The  performance  of  adsorption  systems  is  affected  by 
the  transport  properties  of  the  solutes  in  the  solid  and 
fluid  phases  and  also  by  equilibrium  adsorption  capacity.  The 
presence  of  other  organic  solutes  or  background  organic 
chemicals  may  retard  or  speed  up  the  adsorption  rates  and 
equilibria.  The  particle  diameter,  shape  and  size 
distribution  can  affect  these  parameters. 

1.4  Objectives  of  the  study  :- 


1 )  To  study  from  batch  studies  if  external  mass  transfer 
rate  and  intraparticle  diffussion  rate  varies  with  particle 
size . 

2)  To  assess  the  performance  of  a  fixed  bed  model  to 
predict  break-through  curves  from  column  studies. 

3)  To  develop  a  procedure  to  calculate  a  mean  diameter 
for  a  particle  size  distribution  and  corresponding  mean 
transport  coefficients  for  use  in  batch  and  fixed  bed 
adsorption  models . 


Chapter  2 
Estimation  of  Particle  Size 

There  are  several  industrial  activities  in  which  though 
the  quantities  are  small,  powdered  materials  fulfill  an 
essential  function  in  a  manufacturing  process.  The  useful 
applications  of  powdered  materials  thus  include  many  of  the 
raw  materials  used  in  civil  engineering,  chemical  and  mining 
industries,  pigments  and  countless  other  industries.  Another 
aspect  concerning  materials  in  a  state  of  fine  division,  is 
the  harmful  effect  of  atmospheric  and  industrial  dusts,  mine 
dusts,  smoke  and  grits  from  boilers  and  furnaces. 

Development  of  the  methods  of  sizing  is  the  30b  of  a 
physicist  or  a  chemist,  but  it  is  the  engineer  who  finally 
has  to  adapt  the  results  to  practical  problems.  There  are 
several  sophisticated  instruments  recently  developed  for  size 
analysis.  Precision  methods  of  control  now  adopted  for  modern 
industrial  manufacturing  processes  have  created  a  need  for 
methods  of  particle  size  analysis  to  a  high  degree. 

2.1  Methods  of  sizing 

The  practical  methods  of  size  analysis  thus  comprise: 

1)  Hicroscopial  measurment 

2)  Sedimentation 
3>  Sieving 

When  the  powdered  materials  used  in  industry  extend  over 
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a  range  of  particle  size  that  is  beyond  the  scope  of  analysis 
by  any  one  of  the  above  methods  particle  shape  and  equivalent 
diameters  must  be  taken  into  consideration. 

2.1.1  Hicroscopial  measurment 

The  unique  feature  of  microscopial  measurment  is  that 
the  particles  are  measured  individually,  instead  of  being 
grouped  statiscally  by  some  process  of  classification.  Thus 
the  weight  of  material  examined  is  inevitably  small,  though 
if  there  is  great  variation  in  size,  the  number  of  particles 
involved  to  ensure  a  representative  count  may  be  very  large. 
The  same  difficulties  exist  in  an  even  more  aggravated  form 
with  the  electron  microscope. 

The  lower  limit  of  particle  size  that  can  be  resolved  by 
the  microscope  using  visible  light  is  of  the  order  of  0.2 
microns.  The  lower  limit  that  can  be  resolved  using  ultra- 
violet light  is  0.03  microns  <  Harold  Heywood,  1947  ). 

2.1.2  Sedimentation 

Sedimentation  is  the  most  frequently  used  process  for 
the   size  analysis  of  particles  below  the  sieving  range. 

If  the  particles  are  initially  distributed  uniformly  in 
a  column  of  fluid  at  rest,  after  the  lapse  of  a  certain  time 
there  will  be  a  density  gradient  or  variation  in 
concentration  along  the  height  of  the  column.  Determination 
of  this  density  gradient  after  a  suitable  time  of 
sedimentation  enables  the  proportion  by  weight  of  the  various 


sizes  of  particle  to  be  calculated.  The  variation  in  particle 
concertation  may  be  determined  in  either  of  the  two  following 
ways : 

1)  By  measuring  the  amount  of  material  that  has  settled 
below  the  depth  h,  where  h  is  the  height  of  the  sedimentation 
column  above  the  sampling  level.  This  method  has  been 
designated  "  cumulative  " . 

2)  By  measuring  at  suitable  intervals  of  time  the 
concentration  at  a  particular  depth  h  below  the  surface.  This 
method  ha3  been  designated  "  incremental  ". 

Cumulative  methods  are  rarely  adopted  nowadays  because  the 
curve  relating  mean  concentration  and  time  of  settlement  must 
be  differentiated  graphically  in  order  to  obtain  the  size 
analysis  curve  of  weight  undersize  against  particle  diameter. 
In  the  incremental  method  hydrometers  are  made  use  of  in 
determining  the  concentration  at  a  given  depth  periodically. 
Errors  result  due  to  the  settling  of  particles  on  the 
shoulder  of  the  bulb  and  also  due  to  the  disturbance  of  the 
suspension  when  the  hydrometer  is  constantly  removed  and  re- 
inserted. Moreover  this  process  is  time  consuming. 

2.1.3  Sieving 

Sieving  is  probably  the  easiest  and  certainly  the  most 
popular  method  of  size  analysis  but  it  is  restricted  to 
powders  having  the  greater  proportion  coarser  than  50 
microns.  For  finer  powders  the  method  is  not  generaly  used 
because   of   the  high  cost  of  producing  sieves   with   uniform 
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small  apertures. 

A  sieve  is  an  open  container,  usually  cylindrical  having 
definitely  spaced  and  uniform  openings  in  the  base.  The 
openings  are  square  when  wire  mesh  is  used  and  circular  when 
the  openings  are  formed  by  punching  holes  in  a  metal  plate. 
By  stacking  the  sieves  in  order  of  descending  aperture  size 
and  placing  the  powder  on  the  top  sieve  and  agitating,  the 
powder  gets  classified  into  fractions.  A  closed  pan,  a 
reciever,  is  placed  at  the  bottom  of  the  stack  to  collect  the 
fines  and  a  lid  is  placed  on  top  to  prevent  loss  of  powder. 
Agitation  may  be  either  manual  or  mechanical,  but  mechanical 
is  invariably  chosen.  Results  are  usually  expressed  in  the 
form  of  a  cumulative  undersize  percentage  distribution  in 
terms  of  the  nominal  apertures  of  the  sieve. 

2.1.3.1  Types  of  sieves 

A  variety  of  sieve  aperture  ranges  are  currently  used, 
the  most  popular  being  the  1)  German  Standard,  2)  A. S.T.N, 
standard  Ell-61,  3)  American  Tyler  series  A.F.N.O.R.,  4)  the 
Institute  of  Nining  Standard,  and  6)  the  British  Standard.  In 
all  our  studies  sieving  was  done  using  the  U.  S.  standard 
mesh  size  sieves  which  had  square  openings. 

2.1.3.2  Sieving  errors 

The  apertures  of  a  sieve  may  be  regarded  as  a  series  of 
gauges  which  reject  or  pass  particles  as  they  are  presented 
at  the  aperture.  The  probability  that  a  particle  will  present 
itself  at  an  aperture  depeds  on  the  following  factors: 


11 

1)  The  particle  size  distribution  of  the  powder. 

2)  The  number  of  particles  on  the  sieve. 

3)  The  physical  properties  of  the  particles  (eg.,  surface). 

4)  The  method  of  shaking  the  sieve. 

5)  The  dimension  and  shape  of  the  particle. 

6)  The  geometry  of  the  sieving  surface. 

Wether  or  not  the  particle  will  pass  the  sieve  when  it  is 
presented  at  the  sieving  surface  will  then  depend  upon  its 
dimension  and  the  angle  at  which  it  is  presented. 

The  size  distribution  given  by  a  sieving  operation 
depends  also  on  following  variables: 

1)  Duration  of  sieving. 

2)  Variation  of  sieve  aperture. 

3)  Wear. 

4)  Errors  of  observation  and  experiment. 

5)  Errors  of  sampling. 

6)  Effect  of  different  equipment  and  operation. 

In  all  our  studies  sieving  was  conducted  for  10  minutes 
making  sure  that  all  the  sieves  were  cleaned  well  to  remove 
particles  of  carbon  which  may  have  stuck  in  between  the 
meshes.  Wet  sieving  was  conducted  when  particles  less  than  75 
microns  were  needed  for  isotherm  studies. 

2.2  Particle  shape  and  eqivalent  diameters. 

Many  powdered  materials  used  in  industry  extend  over  a 
range   of  particle  size  that  is  beyond  the  scope  of   analysis 
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by  any  one  method.  for  example  sieving  may  be  necessary  for 
the  coarse  particles,  followed  by  sedimentation  or 
microscopial  measurement  for  subsieve  particles.  The 
equivalent  diameters  corresponding  to  these  various  processes 
are  only  equal  in  the  case  of  spheres.  and  the  relationship 
for  materials  such  as  mineral  particles  will  be  affected  by 
particle  shape. 

Definition  of  particle  shape  as  applied  to  a  number  of 
particles  are  statistical ,  for  in  any  normal  powder  there 
exists  a  wide  variety  of  shapes  ranging  from  flakes  to 
elongated  particles  and  including  some  particles  that  are 
equi -dimensional .  The  numerical  values  of  shape  coefficients 
depend  on  the  characteristic  dimension  chosen  to  represent 
the  particle  size_.J  It  is  also  necessary  to  distinguish 
between  the  effects  of  geometrical  shape,  i.e.  the  closeness 
with  which  the  particle  approaches  the  shape  of  one  of  the 
geometrical  solids,  such  as  sphere,  cube,  tetrahedron,  etc., 
and  the  proportions  of  the  particle,  i.e.,  length,  breadth 
and  thickness. 

Precise  definition  of  the  terms  length,  breadth, 
thickness  is  important.  If  the  particle  is  assumed  to  be 
resting  on  a  plane  in  the  position  of  greatest  stability,  the 
the  breadth  is  the  minimum  distance  between  two  parallel 
lines  tangent  to  the  profile  of  the  particle  when  viewed 
perpendicularly  to  the  above  plane.  The  length  is  defined  as 
the  maximum  distance  between  two  parallel  lines  tangent  to 
the   outline  as  defined  above  and  perpendicular  to  the   lines 
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defining  the  breadth.   Thickness  is  the  distance  between   the 
two   planes   parallel   to  the  plane    defining   the   greatest 
stability  and  tangent  to  the  surface  of  the   particle 
(  Allen,  1974  ). 

If  the  length,  breadth  and  thickness  are  represented  by 
L,  B  and  T,  then  the  flakiness  and  elongation  of  the  particle 
can  be  defined  by  the  following  ratios  : 

Flakiness,  m  =  B/T 
Elongation,  n  =L/B 

If  we  consider  a  particle  to  be  circumscribed  by  a 
rectangular  parallelpiped  of  dimensions  L*B»T,  then  the 
projected  area  of  the  particle.  A,  can  be  represented  by 

A  =  (  PI/4  )*dp2  ■  o^BL (1) 

where    c<a    is    the    area    ratio. 

The    volume    V    of    the    particle    can    be   representd    by 

V    =    c<v,adp3    =    otjBLpj-T (2) 

where 

Pr  is  the  prismoidal  ratio  and  is  equal  to  T/3 

and  °fya  is  volume  shape  coefficient. 

Heywood   (   1947   )   classified   particles   into   tetrahedal , 

prismoidal,  subangular,  and  rounded.  From  his  studies  he  came 

up   with   values  of  04,  and  pr  which  are  shown  in   table   2.1. 

Combining  equations  1  and  2  the  volume  shape  coefficient 
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o(v,a  can  be  directly  given  by  the  following  equation 

*  =        ^        ^_ 

w      A. 

When  the  particle  becomes  equidimensional,  ie.  when  L=B=T  and 
n  =  m  =  l  then  the  volume  shape  coefficient  of  the 
equidimeneional  particle  is  given  by 

^c       =■     ■    — — 


where    0(9    ia    the   volume   shape    coefficient    of    the 

equidimensional  particle. 

Hence 

V,*. 

Hence  we  see  that  shape  can  be  given  as  a  combination  of  the 
relative  proportions  L,  B,  and  T  and  the  geometrical  form  it 
falls  into.  (As  a  result  of  this  we  see  that  there  exists  a 
method  of  seperating  the  effects  of  geometrical  shape  and 
relative  proportions  of  the  particle.' 

Surface  of  a  particle  can  be  given  by 

S  =  <*s,adp2 

where  tfs , ^  is  the  surface  shape  coefficient  of  the  particle. 

Heywood    developed    a   relationship    for    °<s  ■ -1    from 
experiments  which  is  given  by 
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«.*«  ^<%-)4/ifeO 


.5, a. 

where   C  depends  on  the  geometrical  form  of  the  particle   and 

is  a  constant.  Table  2.2   shows  the  values  of  C   and  °6g,a  for 

various    geometrical   shapes   and  irregular   forms   of   the 
particles. 


2.3   Particle  size  distribution 

The  granular  activated  carbon  was  sieved  through  a 
series  of  U.S.  standard  sieves  ranging  from  sieve  no.  12  to 
sieve  no.  40.  Eight  different  size  fractions  where  collected. 
The  particle  size  distribution  of  the  carbon  supplied  by  the 
manufacturer  based  on  weight  is  shown  in  Table  2.3. 

Examination  of  the  particles  using  optical  scanning 
techniques  indicate  these  particles  are  generally  not 
spherical.  (  Particle  size  and  shape  have  a  significant  impact 
on  the  estimated  transport  coefficients  used  in  adsorption 
models. I  Particle  size  and  shape  modification  have  also  been 
shown  to  improve  performances  in  fixed  and  fluidized  beds 
using  molecular  sieve  adsorbents  and  catalyst  particles 
<  Moharir  et .  al . ,  1980  ). 

Activated  carbon  samples  were  examined  using  a  Quantimet 
720-23A-  programmable  and  computerized  image  analyser.  This 
equipment  uses  a  high  resolution  low  noise  video  scanning 
system,   and   is  connected  to  a  PDP-11  minicomputer  for   data 


Table  2.1  Values  for  at   and  pr  for  various  shapes 


Table  2.2  Values  of  °ltsa    and  C  for  various  geometrical 
forms  and  also  for  irregular  particles. 


Shape  group  «oa 


4.36 
2.55 
1.86 


3.3 
3.0 


2.6 

2.1 


Geometrical  forms 

Tetrahedral 

0 

.328 

cubical 

0 

.696 

spherical 

0, 

.524 

Approximate  forms 

Tetrahedral 

0, 

.38 

Angular 

Prismoidal 

0 

.47 

subangular 

0 

.51 

rounded 

0, 

,54 
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Shape  group 

«a 

Pr 

Tetrahedral 

0 

.5  -  0.8 

0 

.4  -  0.53 

Angular 

Prismoidal 

0 

.5  -  0.9 

0 

.53  -  0.9 

Sub-angular 

0. 

.65  -  0.85 

0. 

.55  -  0.8 

Rounded 

0, 

.72  -  0.82 

0 

.62  -  0.75 
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Table  2.3  Particle  size  distribution  in  the 
given  carbon. 


Sieve  size  Percent  fraction 

#12-#14  9.4 

#14-#16  21.0 

#16-#18  17.6 

#18-#20  20.0 

#20-#25  14.5 

#25-#30  9.8 

#30-#35  4.6 

remaining  3 . 1 
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processing . /Parameters  measured  using  this  instrument  include 
particle  area,  perimeter  and  aspect  ratio.n  jThese  parameters 
for  individual  particles,  were  combined  for  a  size  fraction 
to  generate  mean  particle  diametersi 

Quantimet  image  analysis  is  a  recognized  tool  for 
direct  and  quantitative  assessment  of  macro  and  micro 
objects.  Quantimet-720  system  was  used  to  perform  rapid  and 
reliable  analysis  of  coal  particle  samples  categorized 
individually  according  to  their  projected  area,  equivalent 
diameters,  perimeter,  circularity,  shape  factor  and  aspect 
ratio. I 

Quantimet-720  utilizes  a  cathode  ray  tube  vidicon.  The 
optical  images  results  in  the  formation  of  a  matching  image 
made  up  of  electrical  charges  varying  with  image  brightness. 
The  charged  image  is  scanned  by  an  electron  beam  to  produce 
an  output  video  signal. 

Conditions  under  which  the  images  were  digitized  had  to 
be  carefully  controlled.  Electronic  digitizing  equipment  have 
to  be  set  to  establish  optimal  conditions  for  detection.  The 
video  image  of  the  sample  is  transmitted  and  displayed  on  the 
TV  monitor.  The  several  modules  process  the  digitized 
information.  Assuming  proper  detection,  the  digitized  image 
or  the  detected  image  is  a  good  approximation  for  the  shape 
and  area  of  the  object. 

I  The    output    of   measurements   is   in   the    form    of 
distribution  graphs  of  above   mentioned  shape  factors.   Mean, 


19 

median,    standard    deviations   of   the   given   group    were 

automatically  calculated. 

-J 

Errors  could  be  caused  by  particles  touching  one  another 
and  this  was  noticed,  especially  in  small  particles.  To 
remove  this  error  those  particles  which  stuck  to  one  another 
were  not  measured.  In  each  field  a  number  of  adsorbent 
particles  were  grouped  to  produce  an  image  that  filled  the 
size  of  the  screen.  The  results  of  Quantimet-720  measurements 
were  plotted  on  log-probability  paper  as  percent  cumulative 
by  weight  less  than  the  stated  size  against  diameter.  I  The 
particle  size  distribution  of  the  sieve  sizes  in  the  range  of 
#12-#14  to  #30-#35  are  shown  in  the  appendix..  The  particle 
size  distribution  shown  in  the  appendix  is  in  the  form  of 
log-normal  distribution.  From  these  plots  it  is  seen  that  for 
any  particular  size  fraction  there  are  particles  with  varying 
diameters.  Table  2.4  gives  the  mean  area  of  the  particle, 
mean  perimeter  of  the  particle,  mean  bulkiness  and  mean 
aspect  ratio  of  the  particles. 

From  the  table  we  see  that  the  mean  aspect  ratio  is 
almost  equal  to  the  mean  bulkiness  of  the  particle.  It  was 
also  seen  that  the  mean  aspect  ratio  of  the  particles  for 
each  size  fraction  almost  remained  a  constant  at  1.3. 

2.4  Determination  of  particle  diameters 

The  average  diameter  of  the  particles  was   determined  in 
three  different  ways.  They  are  as  follows  :- 
1)  Geometric  mean  diameter  (  Dgm  ). 
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Table  2.4  Mean  properties  of  particles  for  individual 
size  fraction  based  on  Quantimet  analysis. 


Sieve 
size 

Mean  area 
of  particle 

Mean  perimeter 
of  particle 

Mean 
bulkiness 

Mean 

aspect 

ratio 

Number 

of 

particles 

(  mm  ) 

(  mm 

) 

(  per   ) 
4*pi*Area 

len/br 

//12-#14 

2.27 

6.07 

1.29 

1.3 

109 

#14-#16 

1.7 

5.25 

1.29 

1.25 

229 

#16-#18 

1.01 

3.9 

1.20 

1.29 

321 

#18-#20 

0.824 

3.5 

1.18 

1.25 

284 

#20-#25 

0.622 

3.18 

1.29 

1.3 

281 

#25-#30 

0.45 

2.6 

1.20 

1.26 

378 

«0-#35 

0.28 

2.05 

1.19 

1.28 

439 
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2)  Length  mean  diameter     (  Dim  ) . 

3)  Sauter  mean  diameter     <  Dsm  ) . 

In  all  our  studies  a  known  weight  of  carbon  was  used. 
Hence  the  diameters  had  to  be  evaluated  based  on  mass.  The 
following  paragraphs  discusses  the  different  diameters  and 
the  process  of  converting  diameters  represented  by  number 
fractions  into  weight  fractions. 

1)  Geometric  mean  diameter  :-  The  geometric  mean  diameter  of 
the  carbon  particle  was  found  by  finding  the  geometric  mean 
of  the  upper  and  lower  sieves  ie. 


Dgm  ■  <  di*d2  > 


0.5 


2)  Length  mean  diameter  :-  The  length  mean  diameter  is 
represented  in  terms  of  number  distribution  by  the  following 
equation  :  <  Foust  et  al.,  1956  ) 


l- ; i*> 


but 


AH.  =  $»0/0 ® 
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The  left  hand  aide  la  the  differential  change  in  sample  mass 
for  a  differential  change  in  particle  dimensions  for  discrete 
intervals. 


A^     m    (fXJt^ @ 

where   the  left  hand  side  is  the  number  of  particles  in   that 
dimension  class. 

*WI-    7# © 

where    the   dinominator   term   ia   the   mass   of   individual 
particles  and  fpy  is  the  shape  factor  for  volume. 

Hence   the  length  mean  diameter  in  terms  of  mass   reduces 

to 

£   e  %  v3  w  £■ 

where   W  is  the  total  mass  of  carbon  taken.   The   dinominator 
reduces  to  one  and  so  the  equation  reduces  to 

k 


•  i 
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where  Axi  is  the  weight  fraction. 

3)      Sauter  mean  diameter     :-  The  sauter  mean  diameter  is 

otherwise  called  the  volume  surface  mean  and  can  be 
represented  in  terms  of  number  distribution  by  the  following 
equation  : 


using   equations  A,   B,   and  C  from  the  previous  page  we   can 
represent  this  equation  in  terms  of  mass 


Z-  D, 


lm 


*<- 


zk 


Table  2.5  gives  the  formulas  for  the  conversion  from 
number  to  mass  or  vice  versa.  Tables  2.6  and  2.7  give  the 
sieve  size  fraction  and  their  different  diameters  based  on 
quantimet  measurments  with  respect  to  area  of  the  particle 
and  perimeter  of  the  particle  respectively. 

From  the  tables  we  see  that  the  Sauter  mean  diameters 
based  on  perimeter  of  the  particle  had  the  largest  value. 
Sauter  mean  diameters  for  both  area  and  perimeter  of  the 
particle  have  largest  diameters  in  their  respective  classes. 
The  Sauter  mean  diameter  is  a  function  of  the  volume  and  the 
surface  area  of  the  particle.  Moreover  in  any  particular  size 
range,  particle  shapes  vary  significantly.  This  results  in 
particles  with  various  linear  dimensions.  The  presence  of 
irregularly  shaped  particles  that  are  extremely  long  or  flat 
may  result  in  particles  with  sizes  larger  than  the  upper 
sieve  size  and  smaller  than  the  lower  sieve  size.  In 
comparison  to  the  Sauter  mean  diameter  the  geometric  mean 
diameter  is  just  a  function  of  the  dimensions  of  the  mesh 
size.  This  could  possibly  be  a  reason  why  the  Sauter  mean 
diameters  are  much  larger  than  the  geometric  mean  diameters. 
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Table  2. £  Carbon  size  fraction  and  diameters  using  area 
of  the  particle. 


Size 

Geometric 

Length 

Sauter 

fraction 

mean  diameter 

mean  diameter 

mean  diameter 

microns 

microns 

microns 

#12-#14 

1540.0 

1710.0 

1752.0 

#14-#16 

1295.0 

1504.0 

1556.0 

#16-#18 

1091.0  . 

1122.0 

1157.0 

#18-«0 

917.0 

1045.0 

1075.0 

#20-#25 

772.0 

888.0 

923.0 

#25-#30 

647.0 

769.0 

782.0 

#30-«5 

543.0 

638.0 

645.0 

2? 


Table  2.7  Carbon  size  fractions  and  diameters  using 
perimeter  of  the  particle. 


Sieve 

Geometric 

Length 

Sauter 

size 

mean  diameter 

mean  diameter 

mean  diameter 

microns 

microns 

microns 

m-#i4 

1540.0 

1920.0 

1965.0 

#14-#16 

1295.0 

1665.0 

1710.0 

#16-#18 

1091.0 

1255.0 

1297.0 

#18-#20 

917.0 

1148.0 

1184.0 

#20-#25 

772.0 

1010.0 

1067.0 

#25-#30 

647.0 

840.0 

867.0 

#30-#35 

543.0 

655.0 

686.0 

Chapter  3 
Review  of  Related  Literature 

Now  that  the  particle  diameters  in  terms  of  weight  have  been 
estimated  it  is  necessary  to  see  which  of  the  diameters  could 
predict  batch  and  column  adsorption  best. 

3.1.  Isotherms 

An  adsorption  isotherm  is  an  expression  of  the 
equilibrium  at  constant  temperature,  between  the 
concentration  of  a  species  on  the  adsorbent  surface  and  the 
concentration  in  bulk  solution.  In  simpler  terms,  isotherm 
indicates  the  maximum  amount  of  solute  a  particular  type  of 
adsorbent  can  adsorb.  Isotherms  can  be  classified  as  to 
whether  they  describe  the  adsorption  of  a  single  solute  or 
mixture  of  solutes  and  if  it  is  monolayer  adsorption  or 
multilayer  adsorption. 

3,1.1.  Single  solute  monolayer  adsorption 

The  two  simplest  isotherm  equations  for  describing  the 
adsorption  of  a  single  species  are  the  irreversible 
(rectangular)  isotherm  and  the  linear  i3otherm  and  they  are 
represented  by  the  following  equations  : 

q  =  Kr  rectangular  isotherm 

q  ■  KiC linear  isotherm 

where    q   is   the   surface   concentration   and   C    is    the 
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concentration  of  solute  in  the  liquid  phase. 

Unfortunately  neither  equation  is  generally  adequate  in 
characterizing  the  adsorption  of  organics  from  aqueous 
solution  onto  activated  carbon.  The  equilibrium  data  of  a 
very  strongly  adsorbing  species  may  in  some  cases  be 
approximated  by  the  irreversible  isotherm.  The  linear 
isotherm  is  most  likely  to  be  valid  at  low  concertations. 

Two  Vv  linear  models,  the  Freundlich  isotherm  and  the 
Langmuir  isotherm  have  been  frequently  employed  to  describe 
equilibrium.   They  are  represented  by  the  following  equations 

q  ■  KCn  Freundlich  isotherm 

q  ■  QbC/ <  1  *  bC  )  --  Langmuir  isotherm 

The  Freundlich  isotherm  was  originally  developed  empirically. 
The  Langmuir  isotherm  was  theoretically  derived  on  the 
assumptions  of  a  monolayer  as  maximum  adsorption,  a  constant 
energy  of  adsorption,  and  no  migration  of  solute  in  the  plane 
of  the  surface.  It  predicts  a  saturation  surface 
concentration,  q  ,  for  high  solution  concentrations  and  it 
reduces  to  the  linear  form  at  low  concentrations. 

In  case  of  adsorption  by  activated  carbon  in  aqueous 
systems,  neither  Freundlich  nor  Langmuir  equations  may 
describe  the  data  satisfactorily  over  a  range  of 
concentrations.  Therefore,  an  empirical  equation  with  three 
parameters  has  been  proposed  <  Redlich  and  Peterson,  1959: 
Radke  and  Prausnitz,  1972;  and  Mathews  and  Weber,  1977  > .  The 
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equation  is  of  the  form 

q  =  AC/C  1*  BC/3  )  &  N<  1 
At  high  concentrations  this  equation  becomes  Freundlich  and 
at  low  concentrations  it  becomes  linear  and  when  #  -  1  it 
becomes  Langmuir  isotherm.  The  three  parameters  A,  3,  &  can 
be  determined  by  statistical  fit  of  the  equation  to  the 
experimental  data . 

Another  three  parameter  model  recently  presented  by 
Jossens  et  al .  (1978)  is 

C  =  <  q/H  )  exp<  Fqp  ) 

or 

In  C  =  In  q  -  In  H  +  #qep 

This  equation  expresses  solution  concentration  as  a  function 
of  surface  concentration,  p  is  a  constant  such  that  O  <  p  <  1 
and  is  related  to  the  distribution  of  energy  sites  on  the 
surface.  H  and  F  are  functions  of  temperature  only  and  are 
constants . 

Another  type  of  isotherm  expressing  solution 
concentration  as  a  function  of  surface  concentration  is  one 
suggested  by  VanVliet  and  Weber  <  1979,1980  ).  The  equation 
is  as  follows: 


wOv****0*) 


c  =  «-. 

•if   32,   03,   «4  are  coefficients  determined  by  a   nonlinear 
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curve   fitting   method   to   minimize  the  sum   of   squares   of 
residuals  between  experimental  and  predicted  values. 

Sweeny  et  al .  <  1982  >,  suggested  a  different  approach 
for  evaluating  adsorption  isotherm  parameters.  It  assumes  the 
amount  of  adsorbent  as  the  independent  variable  and  final 
concentration  as  the  dependent  variable.  The  equation  is 
given  as  follows: 

c2*(l/b-Co-Cn*q0X)C+CoCn-Co/b-Q°CnX  =  0 

For  Langmuir  isotherm 

(C-CoMC-Cn)-1/n  ♦  KfX  =  0 

For  Freundlich  isotherm 
where 

b   =  parameter  in  the  Langmuir  equation 
C   =  concentration  of  total  organic  carbon 
Cn   ■  total  organic  carbon  resistant  to  adsorption 
Crj    ■  initial  concentration  of  total  organic  carbon   during 

isotherm  tests 
kf   =  parameter  in  the  Freundlich  equation 
n   =  parameter  in  the  Freundlich  equation 
X   ■  concentration  of  activated  carbon 

In  our  studies  isotherm  experimental  data  was  fitted  by 
the  three  parameter  isotherm  equation  as  proposed  by  Mathews 
and  Weber. <1975)  . 

3.1.2  Single  solute  multilayer  adsorption 

The   Brunauer.   Emmet,   Teller   (BET)   model   represents 
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isotherms  reflecting  multilayer  adsorption.  The  BET  model 
assumes  that  a  number  of  layers  of  adsorbate  molecules  form 
at  the  surface  and  that  the  Langmuir  equation  applies  to  each 
layer.  A  further  assumption  of  the  BET  model  is  that  a  given 
layer  need  not  be  completed  prior  to  the  initiation  of 
subsequent  layers.  The  BET  isotherm  reduces  to  the  Langmuir 
model  when  the  limit  of  adsorption  is  a  monolayer.  The  BET 
model  equation  is  represented  as  follows: 

qe=  BCQ°/(<CS-C) <1*CB-1) <C/CS) ) 

where  Cs  is  the  saturation  concentration  of  solute.  C  is  the 
measured  concentration  in  solution  at  equilibrium.  0°  is  the 
number  of  moles  of  solute  at  the  adsorbed  phase  per  unit 
weight  of  adsorbent  in  forming  a  complete  monolayer  on  the 
surface.  qe  is  the  number  of  moles  of  solute  adsobed  per  unit 
weight  at  concentration  C,  and  B  is  the  constant  of  the 
energy  of  interaction  with  the  surface.  This  equation  can 
also  be  put  in  the  linear  form  as  follows  : 

C/<  <Cs-C)qQ)=  1/(B0°)*(B-1/<BQ°)) (C/Cg) 

3.1.3  Multisolute  equilibrium  models 

The   Langmuir   model   of  competitive   adsorption   is   as 

follows : 

qi=<  <0ibiCi)/<l*2b:Cj)) 
J--» 

This   was   first   derived   by   Butler   and   Ockrent   and   the 
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assumptions   were   based   on  the  Langmuir   model   for   single 

solute  systems.   The  surface  concentration  qi  of  adsorbate   i 

is   expressed   as  a  function  of  the  solute  concentrations   of 

all  k  species  in  the  mixture.   The  constants  in  the   equation 

are  obtained  from  single  solute  systems. 

Schay   et  al .   (1957)  modified  the  Langmuir   competitive 

equation  to  improve  its  description  of  equilibrium  data.   The 

model  is  as  follows  : 

K 
qi  =  Qibi<Ci/f7i)/(l*S  bjCa/Tjj) 
>' 

The   interaction   terms  t)j    are  evaluated  by   correlating   the 

model  to  multi-solute  experimental  data. 

Jain   and   Snoeynik  (1973)  introduced  a   bisolute   model 

which   is  a  modification  of  the  Langmuir  and  is  to  be  applied 

to  those  systems  where  Qi  is  not  equal  to  Q2-  The  equation  is 

as  follows  : 


°V,  = 


_C°r<Ok.c.  Qab,C, 


%  -  l  +  bn^+b^Cj. 

It  is  a  semicompetitive  model  where  the  quantity  of  sites 
corresponding  to  Qi  -  0,2  <  where  Qx  >  °<2  *  on  assumed  to  be 
due  to  either  chemical  specificity  or  molecular  exclusions, 
receptive  only  to  solute  1  whereas  the  solutes  compete  for 
the  remaining  sites  corresponding  to  92- 

Mathews   (1975)  extended  the  three-parameter  isotherm  to 


mui tiaolute  adsorption  by  forwarding  the  following  equation  : 

The   constants  are  evaluated  from  single  solute  systems   that 

are   described  by  the  three-parometer  isotherm.  He   proposed 

another   equation   by  applying  the  correlation  procedure   of 
Schay  et  al .  (1957) 


i  +  l^C^V 


Fritz      and   Schlunder   (1974)  have   proposed   a   general 
empirical  multisolute  model  : 

"CO  *-£ 


■Vl- 


fc  -  bcj 


Ei*i«McJ 


For  certain  values  of  the  constants  ,  the  model  reduces  to 
the  Langmuir  competitive  model  or  Mathew's  multisolute 
extension  of  the  three-parameter  isotherm.  Some  of  the  model 
constants  need  to  be  obtained  from  multisolute  data. 

The  ideal  adsorbed  solution  (IAS)  model  has  been  derived 
from  thermodynamics  by  Radke  and  Prausnitz  (1972).  It 
predicts    multisolute    equilibrium    using     single-solute 
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Isotherms   and   is   not  restricted  to  any   specific   type   of 
isotherm.  The  assumptions  in  the  model  are  : 

1)  The  adsorbent  is  thermodynamically  inert 

2)  The  available  surface  area  is  identical  for  all  solutes 

3)  The  liquid  soluton  is  dilute 

4)  The  adsorbed  phase  forms  an  ideal  solution 
The  working  equations  are  as  follows  : 


qT    =    qi    +    q2  (1) 


%■  ■  f  2L  tc,- .  J*  JL  <ti 


1/qT  *  2i/qi°  * 

( 

1  ■ 

-  Zl  )/q2° 

Cl°  =  f(qi°) 

<5) 

C2°  ■  f<q2°> 

(6) 

Ci       ■  Ci°zi 

<7> 

C2        =  C2°C1- 

■*i 

) 

(83 

X 


C4) 


•<D 


where 

qx°  is  the  pure  component  adsorption  for  solute  one,  and 

q2°  is  the  pure  component  adsorption  for  solute  two 

qi°  =  aCi°/<  1  *  BCi   >  (9) 

q2°  ■  aC2°/<  1  *  BC2   >  <10) 
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The  values  of  a ,  B,  and  /g  are  determined  from  single  aolute 
isotherm  studies.  The  values  for  Ci°  and  C20  con  be  solved 
simultaenously  by  substituting  the  equations  9  and  10  into 
the  equations  3  and  4.  From  equations  7  and  8  the  solution 
concentration  Ci  and  C2  can  be  calculated  and  using  the 
values  for  Ci  and  C2»  qi  and  q2  can  be  determined. 

There  are  several  possible  drawbacks  to  the  IAS  model. 
Single-solute  data  are  required  at  very  low  concentrations 
for  the  accurate  calculation  of  spreading  pressure. 
Considerable  computational  effort  is  required  to  solve  the 
simultaneous  equations,  especially  if  the  spreading  pressure 
integration  C  equation  3  )  cannot  be  accomplished 
analytically.  Finally,  the  assumptions  of  an  ideal  adsorbed 
phase  may  fail  at  high  surface  loadings,  leading  to  poor 
predictions. 


3.2  Adsorption  kinetics 

Adsorption  kinetics  involve  determining  the  rate  of 
adsorption  on  the  adsorbent.  The  rate  of  adsorption  is 
determined  by  one  or  more  diffusional  steps.  The  steps  are  as 
follows  : 

1)  Mass  transfer  from  solution  phase  to  the  external  surfaces 
of  the  adsorbent  through  the  liquid  film  surrounding  the 
particle  and  is  known  as  film  transfer;  2)  Transport  of  the 
liquid  within  the  particle  and  is  known  as  intraparticle 
diffuion:  3)  adsorption  reaction  at  the  surface. 
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The  intraparticle  diffusion  can  be  further  divided  into 
pore  and  surface  diffussion.  Pore  diffussion  is  when  the 
solute  is  transferred  to  the  centre  of  the  particle  by  the 
liquid  within  the  pores  and  surface  diffusion  is  when  the 
solute  is  transferred  to  the  centre  of  the  particle  by  moving 
along  the  surface  of  the  adsorbent  and  can  be  said  to  be  a 
result  of  the  concentration  gradient  within  the  individual 
adsorbent  particle.  The  slowest  of  the  three  steps  is  the 
rate  limiting  step.  The  reaction  step  is  usually  very  fast 
and  is  neglected  as  the  rate  limiting  step. 

3.2.1  Single  solute  batch  kinetics  modeling 

Several   models  have  been  developed  taking  into   account 
one   or   more  different  rate  processes.  There  have  also   been 
attempts   to   simplify   these  rate  processes  by   taking   into 
consideration  different  assumptions.  They  are  as  follows  : 
1)   An   overall   or   effective   reaction   rate   expression 
(Thomas,   1944;   Hiester   and  Vermeulen,   1952;   Keinath   and 
Weber,   1968   )  .   2)   A   linear  or   quadratic   driving   force 
approximation   for   intraparticle   diffusion   (Glueckauf   and 
Coates,   1947;   Vermeulen,  1953;  Vassiliou  and  Dranoff,  1962; 
Hall  et  al.  1966;  Cooney  and  Strusi,  1972;  Hsieh  et  al .  1977) 
3)   Film   transfer   as   the  only  rate   controlling   factor 
(Gariepy  and  Zwiebel,  1971;   Zwiebel  et  al .  1972;  Kyte,  1973; 
Keinath,  1977  ) 

Weber    and    Rumer    (1965)   have   proposed    a    model 
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considering  pore  diffusion,  cylindrical  particle  shape  and 
Langmuir  isotherm  for  adsorption  of  benzenesulf onates  on 
actvated  carbon.  Snoeyink  and  Weber  <1968)  proposed  a  similar 
model  but  having  spherical  particles. 

Suzuki  and  Kawazoe  (1974)  proposed  a  model  which  considered 
surface  diffusion  as  the  rate  controlling  step.  The 
dif f usivities  were  found  for  fifteen  different  organics  and 
correlated  to  the  ratio  of  the  boiling  point  of  the  solute  to 
the  adsorption  temperatures. 

Neretnieks  (1967)  has  presented  and  solved  equations  for 
pore  surface  and  combined  diffusion  along  with  film  transfer 
considerations.  The  isotherms  used  were  of  Freundlich  and 
Langmuir  forms. 

Digiano  and  Weber  (1973)  have  considered  film 
resistance,  pore  diffusion  and  Langmuir  isotherm  in  their 
infinite  batch  model  to  describe  data  for  p-nitrophenol  and 
2.4-dinitrophenol  with  activated  carbon.  Infinite  batch 
system  is  one  in  which  a  boundary  condition  of  constant 
concentration  of  solute  in  solution  is  implied.  It  was  found 
that  the  anionic  forms  of  the  solute  did  not  diffuse  as 
rapidly  as  in  their  neutral  forms. 

Hathews  and  Weber  (1977)  developed  the  homogenous 
solid  phase  diffusion  model  with  the  three  parameter  isotherm 
and  succesfully  predicted  adsorption  for  four  solutes 
exhibiting  widely  different  adsorption  equilibrium  and  rate 
characteristics . 
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3.2.2  Single  aolute  kinetics  for  fixed  beds 

Several  models  have  been  developed  using  the  rate 
process  described  earlier  to  predict  breakthrough  profiles  of 
fixed  beds. 

Morton  and  Murril  (1967),  Stuart  and  Camp  (1973)  have 
used  fixed  bed  models  accounting  for  film  and  surface 
diffusion  and  a  nonlinear  isotherm.  The  experimental  data 
fitted  well  to  the  model.  A  similar  model  but  with  axial 
dispersion  was  used  by  Colwell  and  Dranoff  (1969)  and  it  was 
seen  that  the  experimental  data  fitted  the  model  well. 

In  1970  Wheeler  and  Middleman  proposed  a  model  which 
incorporated  three  types  of  resistances.  Intraparticle 
transport  was  described  by  surface  diffusion  and  the  surface 
reaction  was  assumed  to  obey  Michaelis-Menton  kinetics. 
Particle  to  fluid  convective  mass  transfer  was  also 
considered . 

Peel  and  Benedeck  (1980)  have  developed  a  model  based  on 
the  knowledge  of  the  internal  structure  of  the  activated 
carbon.  Here  the  carbon  particle  is  divided  into  two  regions, 
micropores  and  macropores.  In  the  macropores  relatively  rapid 
diffusion  and  adsorption  takes  place  and  the  rest  of  it  takes 
place  in  the  micropores  till  equilibrium  is  reached.  This 
model,  however  did  not  show  any  significant  change  in  the 
model  predictions. 

Weber  and  Pirbazari  (  1982  )  studied  the  adsorption 
characteristics   of   benzene,    p-dichlorobensene,     carbon 
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tetrachloride,  dieldrin,  and  two  PCB's  in  water.  They  used 
the  Michigan  Adsorption  Design  and  Application  Model  (MADAM) 
for  data  synthesis  to  simulate  and  predict  how  carbon  beds 
perform  in  removing  compounds  under  water  treatment 
conditions. 

3.2.3  Single  solute  fixed  bed  models  with  two  or  more 
adsorbents. 

Many  researchers  have  tried  to  use  different  adsorbents 
rather  than  one  adsorbent.  Smith  et  al .  (1959)  studied  the 
adsorption  of  2, 4-dichlorophenol  with  four  different  types  of 
activated  carbon.  They  tried  to  correlate  kinetics  to  the 
carbon  charachteristics  like,  pore  size  distribution  and 
surface  area.  A  relationship  was  apparaent  between  adsorption 
rate  and  micropores  which  were  less  than  250  A  radius. 

Suzuki  and  Kawazoe  (1974)  have  proposed  a  model  for  the 
adsorption  of  2-dodecyl  benzene  sulfonate  on  four  different 
types  of  activated  carbon,  assuming  a  rectangular  isotherm. 

Westermark  (1975)  modeled  four  carbons.  The  water  used 
was  treated  by  sedimentation,  biological  oxidation  and 
coagulation  with  alum.  It  was  further  filtered  to  remove 
suspended  matter.  The  contaminants  were  measured  in  terms  of 
COD.  The  experimental  data  was  corrected  for  nonadsorbable 
and  extremely  adsorbable  compounds  by  varying  the  values  of 
the  isotherm  constants  in  the  linear  isotherms,  and  the  pore 
dif f usivity . 

Holzel  et  al .   (1979)  found  a  fairly  linear  relationship 
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between  the  surface  dif £ usivitiea  for  p-nitrophenol  and  the 
pore  volume  of  the  activated  carbons.  The  p-nitrophenol 
exhibited  a  higher  diffusivity  in  a  carbon  with  larger  pore 
volume . 

Van  Vliet  and  Weber  (1979)  used  two  activated  carbons 
and  eight  synthetic  (carbonaceous  and  polymeric  )  adsorbents 
to  study  the  adsorption  of  p-chlorophenol  and 
p-toluene  sulfonate.  The  adsorbents  differed  widely  in 
capacity,  film  transfer  coefficients  and  surface  diffusivity 
values.  Film  transfer  coefficients  for  synthetic  adsorbents 
could  not  be  estimated  accurately  by  using  literature  values. 
They  therefore  incorporated  particle  shape  and  surface 
topography  factor  in  film  transfer  coefficients. 

Lee  et  al.  (1980),  have  developed  a  model  to  predict 
adsorption  of  humic  substances  on  different  carbons.  The 
model  predicted  correctly  the  fixed  bed  performance  by  using 
film  transfer  coefficients  and  surface  dif f usivities 
independent  of  fixed  bed  experimental  data.  They  also  found 
that  alum  coagulation  before  adsorption  increased  service 
time  of  all  carbons. 

3.2.4  Multisolute  Modelling  for  Batch  Systems 

Mathews  (  1975  )  has  applied  the  homogenous  solid  phase 
diffusion  model  for  single  solute  and  multisolute  adsorotion 
in  a  finite  batch  reactor.  The  three  parameter  isotherm  of 
Radke   and  Prausnitz  was  extended  for  a  bisolute  system  using 
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bisolute  equilibrium  data.  Film  transfer  coefficients  and 
surface  dif f usivities  were  estimated  by  a  parameter  search 
technique  to  fit  single  solute  rate  data  and  these  values 
were  used  to  predict  bisolute  batch  reactor  data.  Runs  with 
a  mixture  of  phenol  and  p-toulene  sulphonate  were  well 
predicted  but  data  fron  phenol  and  p-bromophenol  were  not  30 
well  predicted.  The  mixture  of  phenol  and  dodecyl  benzene 
sulphonate  was  poorly  predicted  and  was  thought  to  be  due  to 
the  interaction  between  the  solutes.  A  change  in  film 
transfer  coefficients  but  the  same  dif f usivities  gave  better 
predictions. 

Liapis  and  Rippin  (  1978  )  proposed  a  model  with  general 
equations  describing  multisolute  adsorption  in  a  finite  batch 
reactor.  The  equations  included  film  and  intraparticle 
resistance.  Intraparticle  resistance  was  due  to  the  pore, 
surface  or  combined  diffusion.  Numerical  solutions  of 
bisolute  adsorption  with  pore  or  surface  diffusion  were  given 
in  the  case  of  independendtly  diffusing  solutes.  A  solute 
could  be  assumed  to  diffuse  into  the  adsorbent  either 
independent  of  the  concentration  or  with  dependence  on  other 
solutes.  Experimental  data  for  2-butanol  and  t-amyl  alcohol 
adsorbed  on  activated  carbon  were  tested  against  the 
predictions  of  the  model.  The  single  solute  data  was 
described  by  the  Freundlich  isotherm,  while  bisolute  data  was 
described  by  the  five  parameter  equation  of  Fritz  and 
Schlunder  (  1980  ).Film  transfer  coefficients  were  calculated 
from   literature  and  surface  dif f usivities  were  estimated   bv 
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trial  and  error.  Both  models  described  the  adsorption  model 
well  and  they  also  found  that  pore  diffusivity  happened  to 
decrease  with  increasing  starting  concentration. 

Adsorption  of  several  phenolic  compounds  on  activated 
carbon  was  tested  by  Fritz  et.  al (  1980  )  to  develop  their 
model  for  single  solute  and  bisolute  systems.  Freundlich 
isotherms  were  used  for  single  solute  data  but  IAS  theory  or 
Fritz-Schlunder  equations  were  used  for  bisolute  data.  Both 
film  and  intrparticle  transport  were  used  in  this  model.  They 
found  that  for  low  concentrations  of  the  order  of  10"4  M,  a 
simple  film  resistance  model  predicted  the  data  well  and  for 
high  concentrations  intraparticle  transport  had  to  be  taken 
into  consideration.  Deviations  in  the  predictions  from 
experimental  data  were  thought  to  be  due  to  the  diffusional 
interaction  between  solutes. 

3.2.5  Multisolute  Modelling  for  Fixed  Beds 

Crittenden  (1976)  and  Crittenden  and  Weber  (1978)  have 
studied  single  solute  and  bisolute  adsorption  of  the  same 
solutes  and  adsorbent  as  used  by  Mathews  (1975) .  They  used 
homogeneous  solid  diffusion  model  including  film  resistance. 
Influent  concentrations  were  varied  20*6.  The  surface 
diffusivities  were  calculated  from  single  solute  batch 
studies,  while  film  transfer  coefficients  were  estimated  from 
literature.  The  model  predictions  were  good  for  single  solute 
data   but   not   for   bisolute   fixed   bed   experiments.    The 
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discrepancies  observed  were  thought  to  be  due  to  experimental 
data  scatter,  poor  equilibrium  description  or  failure  of  the 
assumption  of  independent  diffusion  of  the  solutes. 

Liapis  and  Pippin  (  1978  )  have  proposed  another  fixed 
bed  model  which  includes  film  resistance,  pore  diffusion  and 
axial  dispersion.  The  model  incorporates  dif f usivities  from 
batch  experiments  and  predictions  for  two  alcohols  were  found 
to  be  in  good  agreement  with  the  experimental  results.  The 
model,  however,  failed  to  describe  the  overshoot  of  the 
weaker ' species.  Also,  increase,  in  bed  length  was  observed  to 
induce  a  larger  displacement  of  this  species. 

Crittenden  et.  al  (  1980  ),  studied  the  adsorption  of  p- 
nitrophenol  and  p-bromophenol  on  a  duolite  A-7  resin  with  a 
model  that  incorporated  homogenous  solid  diffusion  and  film 
transfer.  Film  transfer  coefficients  were  estimated  from 
literature.  The  equilibrium  data  was  described  by  Langmuir 
competitive  isotherm.  The  surface  dif f usivities  were 
calculated  by  fitting  the  model  to  experimental  fixed  bed 
data  for  simultaneous  feeding  of  two  solutes.  Good  agreement 
was  observed  between  predictions  of  the  model  and  the 
experimental  bed  run  data  for  sequential  feeding  of  two 
solutes. 

Famularo  et  al  (  1980  >  ,  have  developed  a  micro/macro 
shell  model  representing  diffusion  through  an  adsorbent 
having  a  bimodal  pore  size  distribution.  An  adsorbent 
particle  was  divided  into  two  regions,  a  spherical  core  with 
micropores   enclosed  by  an  annular  sheil  with   macropores.   A 
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linear  driving  force  was  assumed  to  be  the  driving  force  for 
the  transport  through  each  region.  Film  resistance  was  also 
incorporated.  Kinetic  parameters  were  determined  for  phenol 
and  p-nitrophenol  with  carbon  as  adsorbents,  by  single  solute 
batch  tests  and  isotherms  were  previously  determined. 

Thacker  et .  al  (  1981  ),  have  studied  adsorption  of  3,5- 
dimethyl  phenol  and  Rhodamine  6G  on  four  different  activated 
carbons.  One  carbon  was  pretreated  with  aqueous  chlorine  to 
evaluate  the  effect  of  chlorine  pretreatment  on  adsorption. 
The  film  transfer  and  homogenous  surface  diffusion  were  the 
mass  transfer  steps  incorporated  into  the  model.  The  single 
solute  data  were  fitted  by  the  Freundlich  and  Myers  isotherms 
whereas  the  bisolute  data  was  fitted  by  the  IAS  model.  The 
model  predicted  well  for  single  and  bisolute  bed  studies  of 
adsorption,  desorption  and  sustained  step  changes  in  the 
influent  concentration. 

Yen  and  Singer  (  1984  )  applied  the  IAS  model  with  a 
modified  calculation  to  test  ten  sets  of  binary  and  ternary 
phenolic  mixtures.  The  Langmuir  competitive  model  was  used 
for  comparison.  Here  they  found  that  the  IAS  model  gave 
better  predictions  than  the  Langmuir  competitive  model  in  all 
cases  of  study. 


Chapter  4 
EXPERIMENTAL  MATERIALS  AND  METHODS 


4.1   Materials 

4.1.1  Solutes 

The  solutes  used  in  the  study  were  phenol  and 
parachlorophenol .  The  properties  of  both  solutes  are  shown  in 
the  table  4.1.  Both  the  solutes  were  obtained  in  the 
crystallised  form  from  J.T.Baker  Chemical  Co;  (Phillipsburg , 
NJ  > .  Phenol  has  a  water  solubility  of  82  gm/1  whereas 
parachlorophenol  has  a  solubility  of  2.70  gm/1. 

For  all  experiments  stock  solutions  were  initially 
prepared  and  in  no  experiment  were  the  crystals  directly 
introduced  into  the  reactors.  All  solutions  were  made  up 
using  tap  water  that  was  passed  through  a  three  foot  carbon 
column  to  remove  organics  and  other  adsorbebles  present  in 
the  water.  The  water  passing  through  this  carbon  column  was 
used  for  equilibrium  experiments,  batch  studies  and  column 
experiments. 

4.1.2  Adsorbent 

The  adsorbent  used  was  activated  carbon  <  Carborandum  ) 
supplied  by  CECA  Inc.,  Pryor,  Oklahoma.  The  properties  of  the 
carbon  used  are  listed  in  table  4.2. 
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Table  4.1  Solutes  and  their  proporties. 
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Phenol 


PCP 


Supplier 

J.T.  Baker  Chm. 

J.T.  Baker  Chm, 

Company 

Company 

Reagent  grade 

Baker  TM 

Baker  TM 

Formula 

Molecular  Weight 

94.1 

128.56 

*Ka 

Solubility  in 
water  per  litre 

82gm  /L 

2/^0gm.  /L 

Diffusivity  @ 

25  deg.  C.  in 

2, 
cm  /sec 

1.04*10-5 

9.9*10-6 

1*8 


Table  4.2  Proporties  of  the  adsorbent. 


Manufacturer 
U.S.  Mesh  size 
Raw  material 


CECA  Inc.,  Pryor,  Okhlahoma 

12/40 

Bituminus  coal 


Physical  Proporties 


2 
Surface  area,  m  /gm 

1000-1100 

Apparent  density,  gm  /cc 

0.47 

Particle  density  wetted 

in  water,  gms/cc 

0.60 

Effective  size,  mm 

1.9  or  less 

H9 
4.2  Methods 

4.2.1  Analytical  Methods 

Solution  concentrations  of  phenol  and  parachlorophenol 
were  analyzed  by  ultraviolet  light  absorption 
spectrophotometry  with  a  Baush  and  Lomb  spectronic  710 
spectrophotometer.  Initially,  tests  were  made  to  find  the 
maximum  wavelengths  of  absorption  for  phenol  and 
parachlorophenol.  Once  these  were  established,  several 
solutions  of  different  concentrations  for  phenol  and 
parachlorophenol  were  prepared  and  these  solutions  were  used 
to  obtain  the  calibration  charts  for  each  solute 
independently  at  its  wavelength  of  maximum  absorbance.  Phenol 
was  analysed  at  a  wavelength  of  263  nM  and  parachlorophenol 
at  a  wavelength  of  279  nM . 

In  the  case  of  the  bisolute  system  Friedel's  method 
(1951)  was  made  use  of  to  compute  individual  concentrations. 
This  method  is  based  on  the  assumption  that  absorbance 
measured  in  a  mixture  is  the  sum  of  the  absorbance  of  each 
compound  at  that  wavelength.  Equations  developed  at  each 
wavelength  is  as  follows: 

where 

A268.  A279  ■  Total  absorbance  at  263  nM  and  279  nM 
CpcP.   ^Phenol   *   Concentration  of   p-chlorcphenol   and 
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phenol  respectively, 
^hsnoil'  ^cPl  -  Molar  absorpti vities  of  pheonol  and 
p-chlorophenoi  at  wavelength  2&8  nM 
^>henol2>  ePcP2  ■  Molar  absorpti vities  of  phenol  and 
p-chlorophenol  at  wavelength  279  nM 
From  the  above  set  of  equations  we  see  that  the  only  unknowns 
are  Cpcp  and  Cpnenoi  and  since  there  were  two  equations  these 
concentrations  were  established. 

4.2.2  Analysis  of  particle  diameters. 

This  has  already  been  discussed  in  chapter  two  and 
table  4.3  and  table  4.4  lists  the  the  diameters  of  the 
particles  based  on  the  area  and  perimeter  of  the  particle 
respectively. 

4.3  Experimental  Methods 

4.3.1  Equilibrium  Experiments 

Equilibrium  experiments  were  conducted  to  determine  the 
capacity  of  carbon  at  various  solution  concentrations.  This 
was  determined  by  the  standard  bottle-point  method.  This 
method  consists  of  allowing  a  solution  of  known  initial 
concentration  to  come  into  contact  with  different  amounts  of 
carbon  doses  until  equilibrium  is  reached.  The  carbon 
used  for  these  experiments  were  the  ones  passing  through  75^im 
and  retained  on  the  53/«  sieve.  This  carbon  was  obtained  by 
wet   sieving.   The  carbon  obtained  was  then  piaced   in   glass 
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Table  4.3  Carbon  size  fraction  and  diameters  using  area 
of  the  particle. 


Size 
fraction 

Geometric 
mean  diameter 

Length 

mean  diameter 

Sauter 

mean  diameter 

microns 

microns 

microns 

#12-#14 

1540.0 

1710.0 

1752.0 

#14-116 

1295.0 

1504.0 

1556.0 

#16-#18 

1091.0 

1122.0 

1157.0 

#18-#20 

917.0 

1045.0 

1075.0 

#20-#25 

772.0 

888.0 

923.0 

#25-#30 

647.0 

769.0 

782.0 

#30-#35 

543.0 

638.0 

645.0 

Table  4.4  Carbon  size  fractions  and  diameters  using 
perimeter  of  the  particle. 
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Sieve 
size 


Geometric 
mean  diameter 


Length 

mean  diameter 


Sauter 

mean  diameter 


microns 


microns 


microns 


*12-#14 

1540.0 

1920.0 

1965.0 

#U-#16 

1295.0 

1665.0 

1710.0 

#16-#18 

1091.0 

1255.0 

1297.0 

#18-#20 

917.0 

1148.0 

1184.0 

#20-#25 

772.0 

1010.0 

1067.0 

#25-#30 

647.0 

840.0 

867.0 

//30-//35 

543.0 

655.0 

686.0 
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beakers  and  dried  at  a  temperature  of  120  deg.C.  To  teat  if 
the  carbon  was  totally  free  of  all  moisture.  a  sample  of  the 
carbon  was  set  aside  in  the  heater  and  when  the  weight  of 
this  carbon  did  not  differ  from  two  consecutive  readings 
taken  on  two  consecutive  days  it  was  presumed  that  the  carbon 
was  free  of  moisture.  It  took  five  days  for  the  carbon  to 
reach  this  state. 

Different  amounts  of  carbon  were  then  added  carefully 
to  marked  bottles  with  teflon  seal  caps.  To  each  bottle  100 
ml  of  solution  of  known  concentration  was  added.  The  pH  of 
the  water  used  to  make  up  the  solution  was  7.  This  was  done 
by  initially  filling  up  a  jar  of  4  liter  capacity  with  tap 
water.  One  thousand  mL  of  tap  water  was  then  filled  into  a 
beaker  and  the  pH  of  this  water  was  measured  while  the  water 
was  being  stirred.  The  pH  of  the  water  varied  from  day  to  day 
studies  but  remained  between  8.5  and  8.7.  Depending  on  the  pH 
of  the  water  a  certain  quantity  of  3 .  SN  sulphuric  acid  was 
introduced  into  the  water  in  the  beaker  by  using  a  burette 
until  the  pH  reached  7.  From  this  the  amount  of  sulphuric 
acid  needed  to  lower  the  pH  of  the  water  in  the  jar  was 
established.  This  water  was  then  made  use  of  to  make  up  the 
solutions.  Once  the  solutions  were  introduced  into  the 
bottles,  the  bottles  were  then  tightly  closed  and  sealed  to 
make  sure  that  they  were  air  tight.  The  bottles  were  then 
placed  in  an  incubator  shaker  <  Labline  Instruments  111  > . 
Since  all  experiments  were  to  be  conducted  at  25  deg.C  care 
was  taken  to  see  that  the  temperature  in  the  incubator  was  25 
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deg.C.  The  shaker  was  set  in  such  a  way  that  the  carbon 
particles  were  suspended  and  no  carbon  stuck  to  the  sides  of 
the  bottles.  Three  bottles  of  solution  with  no  carbon  were 
kept  aside  in  the  incubator  shaker  to  see  if  the  glass 
adsorbed  any  of  the  solutes  or  if  there  was  any  change  due  to 
biological  reaction.  Three  bottles  of  solution  of  the  same 
concentration  and  the  same  weight  of  carbon  were  set  aside  to 
see  when  equilibrium  was  reached.  Both  phenol  and 
parachlorophenol  reached  equilibrium  after  three  days. 

Once  equilibrium  was  reached,  the  samples  were  filtered. 
The  filtrate  was  then  analyzed  for  the  remaining 
concentration . 

The  surface  concentration  qe  in  equilibrium  with  the 
solution  concentration  was  then  computed  using  the  following 
set  of  equations 


it  .&-O& 


w 

where, 

qe   ■  surface  concentration  in  mmoles/gm  or  mg/gm 
Co  -    initial  concentration  in  mmoles/L  or  mg/L 
CQ   -  final  concentration  in  mmoles/L  or  mg/L 

V     =  volume  of  the  solution  in  L  <  liters  > 

W     =  weight  of  adsorbent  in  gm 

4.3.2  Batch  Experiments 

The  carbon  for  batch  experiments  was  obtained  by  sieving 
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the  desired  fractions.  Each  size  fraction  was  washed 
thoroughly  with  deionized  water  to  remove  fines.  The  samples 
were  then  dried  in  an  oven  at  120  deg.C  until  constant  weight 
was  attained.  The  carbon  was  then  stored  in  glass  bottles  and 
placed  in  an  incubator. 

Experiments  were  conducted  in  a  rectangular  plexiglass 
vessel  with  internal  dimensions  LxBxH  of  34x30x35  cm.  The 
schematic  setup  of  the  batch  experiment  is  as  shown  in 
Figure  4.1.  For  all  the  experiments  24  litres  of  solution  of 
known  concentration  was  made.  The  water  in  the  vessel  was 
stirred  by  a  four  bladed  impeller.  The  temperature  was 
maintained  at  25  deg.C  from  the  beginning  and  until  the  end 
of  the  experiment.  Before  the  experiment  could  be  run  the  pH 
of  the  water  in  the  vessel  was  maintained  at  7.  This  was  done 
by  using  the  same  approach  used  for  isotherm  studies,  and  by 
making  use  of  3.6N  sulphuric  acid. 

Once  this  was  set  up  a  known  weight  of  carbon  was  added. 
The  speed  of  the  impeller  was  maintained  at  a  speed  of  800 
rpm  for  particles  in  the  range  of  12-14  to  18-20,  and  at  a 
speed  of  700  rpm  for  particles  in  the  range  of  20-25  to  30- 
35.  The  lower  speed  could  not  be  used  for  the  larger  sizes 
because  the  particles  were  not  completely  suspended.  Table  4- 
5  shows  the  size  fraction  and  the  impeller  speed  used  for  the 
different  size  fractions.  The  same  impeller  speeds  were  used 
for  both  phenol  and  parachlorophenol  experiments.  In  the  case 
of  mixture  of  size  fractions  even  if  one  of  the  sizes  fell  in 
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Table  4.5  Impeller  speeds  for  different  size  fractions 
for  batch  studies. 


Sieve  size  Impeller  speeds 

#12-#14  800 

#14-#16  800 

#16-#18  800 

#18-#20  800 

#20-#25  700 

#25-#30  700 

#30-#35  700 
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the  higher  impeller  rpm  range  the  higher  impeller  speed  was 
used  for  the  experiment. 

After  the  addition  of  carbon  readings  were  taken  at 
intervals  of  every  minute  for  ten  minutes  and  then  readings 
were  taken  every  fifteen  or  twenty  minutes  for  four  hours. 
The  sample  was  collected  using  a  clean  pipette  with  a  teflon 
filter  attached  so  that  no  carbon  entered  the  sample.  The 
graphs  of  C/Co  were  then  plotted  against  time. 

4.3.3  Fixed  bed  experiments 

The  experimental  set  up  for  the  fixed  bed  experiments 
are  as  shown  in  the  Figure  4.2.  Two  sizes  of  columns  were 
used.  For  the  range  of  particles  12-14  to  18-20  the  column 
with  a  diameter  of  5.08  cm  was  used,  and  for  the  particles  in 
the  range  of  20-25  to  30-35  the  column  with  a  diameter  of 
3.1  cm  was  used.  Table  4-6  shows  the  operating  conditions  for 
fixed  bed  experiments  for  both  phenol  and  parachlorophenol . 
In  the  case  of  a  mixture  of  sizes  even  if  one  of  the  carbon 
size  fractions  fell  in  the  range  of  the  bigger  column 
diameter  then  the  operating  conditons  for  the  column  with  the 
bigger  diameter  was  used. 

The  carbon  in  both  the  columns  was  supported  using  a 
brass  wire  mesh.  The  carbon  for  fixed  bed  studies  was 
prepared  by  sieving  different  size  fractions  and  washing  each 
size  fraction  thoroughly  in  deionized  water  to  remove  fines. 
This  carbon  was  then  dried  in  an  oven  at  120  deg .  C  for 
sufficient  time  till  constant  weight  was  attained.   The  dried 
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Table  4.6   Fixed  bed  operating  conditions 


Sieve  size     Column  diameter   Flowrate   Hydraulic  loading 
fraction  rate 


#12-#14 

5.08 

fJA-416 

5.08 

#16-#18 

5.08 

#18-#20 

5.08 

#20-123 

3.1 

#25-#30 

3.1 

#30-#35 

3.1 

ml/min       m/min 


500 

0.247 

500 

0.247 

500 

0.247 

500 

0.247 

125 

0.166 

125 

0.166 

125 

0.166 
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samples  were  then  stored  in  glass  jars  in  an  incubator  at 
constant  temperature. 

The  column  was  fed  with  tap  water  from  a  tank  in  which 
constant  head  was  maintained  by  overflow.  The  water  collected 
into  this  tank  initially  passed  through  a  three  foot  carbon 
column.  Two  different  dimeter  of  columns  were  used  depending 
upon  the  size  range  of  the  particle.  For  the  smaller  diameter 
column  two  variable  flow  pumps  were  manipulated  to  force  the 
water  as  well  as  the  concentrated  solution  into  the  column 
at  a  total  flowrate  of  125  mL/min.  The  flowrate  of  the 
solution  was  12.5  mL/min.  For  the  bigger  column  a  variable 
flow  master  flex  pump  was  used  along  with  one  variable  flow 
pump  to  force  the  water  as  well  as  the  solution  into  the 
column  at  a  flow  rate  of  500  mL/min.  The  flowrate  of  the 
solution  in  this  case  was  50  ml/min.  The  solution  was  stored 
in  a  ja_r._gf  20  liter  capacity.  The  solution  was  maintained  in 
the  jar  at  a  concentration  10  times  higher  than  what  was 
required  in  the  experiment  and  the  flow  from  the  jar  was 
adjusted  so  as  to  have  proper  dilution  to  the  desired 
concentration  in  the  system. 

Predetermined  amount  of  carbon  was  introduced  into  the 
column.  In  case  a  mixture  of  sizes  was  used  then  the  largest 
size  was  introduced  first  followed  by  the  next  smaller  size 
and  so  on  so  that  the  bed  was  completely  stratified  according 
to  size.  Water_j»as  then  passed  through  the  column  for_  about 
half  an  hour  and  the  walls  of  the  column  were  gently  tapped 
to  remove  all  the  entrapped  air. 
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Before  the  column  could  be  operated  the  temperature  in 
the  constant  head  tank  was  maintained  at  a  temperature  of 
about  30  deg.  C  so  that  when  this  water  mixed  with  the 
solution  from  the  jar  the  temperature  of  the  water  entering 
the  column  was  exactly  25  deg.  C.  Once  again  the  pJL_  was 
maintained  at_7_.  For  this  the  solution  in  the  jar  was  made  up 
with  -  tap  water  maintained  at  pH  7.  For  all  column  runs  the 
f lowxats  of  water  from  the  tank  was  9  times  higher  than  that 
from  the  solution  jar.  The  pH  of  the  water  in  the  tank  was 
measured  and  the  amount  of  3.6N  sulphuric  acid  required  to 
bring  down  the  pH  to  7  was  also  noted.  Since  the  flowrate 
from  the  tank  was  9  times  higher  than  that  from  the  solution 
jar  ninei  times  more  acid  was  to  be  addded  in  the  solution 
jar  to  maintain  the  pH  of  the  water  entering  the  column  at  7. 
This  procedure  was  adopted  for  maintaining  the  pH  at  7  in  the 
column . 

The  colujnn  was  operated  almost  sixteen  hours  after  the 
initial  setup,  i.e.  until  the  column  was  almost  95« 
saturated.  Effluent  readings  were  taken  at  intervals  of  an 
hour — and  were  analyzed  for  effluent  concentrations.  The 
nondimensional  values  of  C/Cq  were  then  plotted  against  time. 


Chapter  5 
Kinetic  Model 


Film  transfer,  surface  diffusion  and  pore  diffusion  are 
important  rate  controlling  steps  in  adsorption  phenomenon. 
Depending  upon  the  conditions,  one  of  these  or  a  combination 
of  these  will  control  the  rate  of  adsorption.  Pore  and 
surface  diffusion  are  two  parallel  processes.  Besides  these, 
the  equilibrium  relatiohship  is  also  important  and  the 
isotherm  shape  has  been  found  to  influence  the  shape  of  the 
breakthrough   curve  <  Chakravorti  and  Weber,  1974  ). 

5.1   Model  Equations 

Assumptions 

The  following  assumptions  were  made  to  determine  the 
equations  for  a  finite  batch  model  as  well  as  a  fixed  bed 
model . 

1)  The  rate  of  adsorption  is  diffusion  limited,  ie., 
local  equilibrium  prevails  at  the  external  surface  of  the 
adsorbent. 

2)  Liquid  diffusion  resistance  occurs  at  the  external 
surface  and  can  be  described  as  film  transfer. 

3)  The  adsorbent  particle  is  a  homogenous  solid. 
4>   The  particle  is  spherically  shaped. 

5)   Plug  flow  occurs  in  the  fixed  bed. 
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5.1.1  Finite  Batch  Reactor 

A  finite  batch  adsorber  is  composed  of  a  finite  volume 
of  solution  in  which  the  adsorbent  particles  are  rapidly 
agitated.  The  finite  batch  adsorber  model  for  a  single  solute 
consists  of  the  following  set  of  equations: 

1)  The  mass  balance  equation  on  the  adsorbent  particle 
is  given  as 


The  first  term  on  the  left  hand  side  of  the  equation  gives 
the  rate  of  accumulation  of  solute  into  the  fluid  in  the 
pores  and  the  second  term  gives  the  rate  of  accumulation  on 
the  surface  of  the  particle.  The  first  term  on  the  right  hand 
side  of  the  equation  gives  the  rate  of  homogenous  solid 
diffusion  while  the  second  term  is  the  rate  of  pore 
diffusion.  If  homogenous  solid  diffusion  is  considered 
important  and  the  fluid  phase  accumulation  term  is  neglected, 
and  since  surface  concenteration  is  greater  when  compared  to 
concentration  in  pore  fluid  for  small  values  of  ep ,  the 
equation  reduces  to 


This   equation   was   developed  by  Rosen   in   1952.   The   film 
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transfer  can  be  taken  into  account  in  the  diffusion  model  by 
applying  appropriate  boundary  conditions.  In  the  homogenous 
solid  diffusion  this  condition  is  obtained  by  equating  the 
rate  of  mass  transfer  through  the  surface  film  to  the  rate  of 
change  of  average  concentration  of  the  particle  and  is 
expressed  as 


The  initial  and  boundary  conditions  are 


ra  fc>/0j  X=0;     |^_0 


5.1.2  Fixed  Bed  Model 

The  following  equations  apply  for  adsorption  in  a  Single 
solute  fixed  bed  adsorber: 

>-/*  -  &t  *  c**  so  © 
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Equations  1  to  5  give  the  rate  of  accumulation  of  solute  on 
the  surface  of  the  adsorbent  particle  through  homogenous 
solid  diffusion.  Ds  is  the  solid  phase  diffusion  coefficient. 
The  radial  distance  from  the  centre  of  the  adsorbent  particle 
is  r.  The  surface  concentration  is  q.  Equation  2  states  that 
the  initial  condition  that  at  time  t=0  the  surface 
concentration  at  any  radial  distance  inside  the  particle  is 
zero.  Equation  3  implies  that  at  any  time  greater  than  zero, 
the  concentration  gradient  at  the  centre  of  the  particle  is 
zero.  The  radius  of  the  particle  is  R.  Equation  4  gives  the 
mass  balance  equation  on  the  adsorbent  particle.  The  rate  of 
mass  transfer  through  surface  film  is  equal  to  the  rate  at 
which  the  average  concenteration  of  the  particle  changes.  The 
film  transfer  coefficient  is  kf,  Cs  is  the  solution 
concentration  at  the  external  surface  of  the  particle  and  p 
is  the  apparent  density  of  the  particle.  Equation  5  is  the 
isotherm  relation  between  the  surface  concentration  and 
solution  of  the  solute  with  qs  as  the  surface  concentration 
at  the  external  surface  of  the  particle.  Equations  6.  7,  and 
8  apply  to  fixed  bed  conditions.  Equation  6  gives  the  mass 
balance  of  the  solute  at  any  point  in  the  bed.  The  rate  of 
change  in  concentration  of  the  solute  is  equal  to  the  sum  of 
the  convective  change  in  concentration  inside  the  bed  and  the 
rate  of  mass  transfer  through  the  film  to  the  solid  surface. 
Here.   V  is  the  superstitial  velocity  in  the  fixed  bed.  "St-,  is 
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the  porosity  of  the  bed  and  Z  is  the  axial  distance  inside 
the  bed.  Equation  7  is  initial  condition  to  the  bed  at  time 
t=0,  when  solution  concentration  inside  the  bed  is  zero. 
Equation  8  implies  that  at  any  time  greater  than  zero,  the 
solution  concentration  at  the  entrance  to  the  bed  is  equal  to 
the  initial  concentration.  Lb  and  Co  are  length  of  the  bed 
and  initial  solution  concentration,  respectively. 

5.2   Solution  of  the  Model 

Generally,  in  the  case  of  a  nonlinear  isotherm,  the 
differential  equations  in  the  adsorbent  model  are  solved 
using  numerical  methods.  Solutions  can  be  obtained  using 
finite  difference,  orthogonal  collocation,  or  finite  element 
methods.  The  present  model  uses  orthogonal  collection  method 
to  solve  the  differential  equations  numerically. 

In  the  collocation  method  the  solution  of  the 
differential  equations  is  approximated  by  a  trial  function. 
This  trial  function  has  constants  and/or  functions  such  that 
when  this  trial  solution  is  substituted  into  the  differential 
equation  the  residual  is  forced  to  zero  at  the  collocation 
points.  Orthogonal  collocation  is  a  special  case  of 
collocation  method,  in  that,  the  trial  functions  are  a  set  of 
orthogonal  polynomials  and  the  collocation  points  are  the 
roots  of  these  polynomials. 

By  the  use  of  collocation  points  the  partial 
differential  equations  are  reduced  to  first  order  ordinary 
differential  equations.  These  ordinary  differential  equations 
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are  solved  by  a  computer  subroutine  program  developed  by 
Gear  (  Gear,  1976;  Hindmarsh,  1974  ).  This  subroutine 
utilizes  the  Adams  predictor-corrector  method  for  numerical 
integration  of  the  ordinary  differential  equations. 

The  program  was  modified  by  Mathews  and  Kulkarni  (1983) 
and  treated  the  whole  bed  as  a  series  of  small  beds 
containing  one  particular  size  of  carbon. 

As  the  first  bed  becomes  saturated,  effluent  from  this 
bed  becomes  influent  for  the  next  bed  at  respective  times. 
The  series  of  effluent  becomes  influent  to  the  next  bed  until 
the  last  layer,  the  effluent  of  which  is  the  final  effluent 
of  the  whole  bed. 

In  the  original  program  single  solute  equilibrium 
constants  were  determined  by  fitting  a  Freundlich  isotherm 
equation  to  the  equilibrium  data.  The  curved  portion  of  the 
data  near  small  concentrations  was  represented  by  a  linear 
isotherm.  However  when  this  was  done  there  were  some  portions 
on  the  isotherm  which  could  not  be  represented  by  a  linear 
fit.  This  resulted  in  erroneous  predictions.  Hence  the 
program  was  modified  by  Mathews  <  1984  )  so  that  equilibrium 
can  be  represented  by  the  three  parameter  isotherm 
requirements . 

There  are  several  input  data  requirements  for  the 
program.  These  are  1) weight  of  carbon  2) length  of  the  bed 
3)radius  of  the  adsorbent  particle  4)molecular  weight  of  the 
solute    5) influent    concentration    6) isotherm     constants 
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7)collocation   conatanta  3) intrapaticle  diffusion  coefficient 
and  9) external  mass  transfer  coefficient. 

The  intraparticle  diffusion  coefficient  was  established 
from  batch  studies.  The  external  mass  transfer  coefficient 
was  established  using  a  method  suggested  by  Dwiwedi  and 
Upadhyay <1977> .  The  set  of  equations  used  to  determine  the 
external  mass  transfer  coefficient  are  as  follows: 

where 

«     ■  porosity  of  the  bed 

Jd    =  mass  transfer  factor 

Nj{e   =  Reynolds  number  <  DpG/^>  ) 

kf    =  film  transfer  coefficient 

u      =  superficial  fluid  velocity  (  G//»  ) 

NSc   =  Schmidts  number  (  U/(/dD[,)) 
here   Dl  is  the  molecular  diffusivity  and  was  determined  from 
the  equation  of  Wilke  and  Chang  <  1955  >,  which  was  given  by 

where 

T  =  absolute  temperature  in  deg .  K 

M  ■  molecular  weight  of  the  solvent 

I*      -  solution  viscosity  in  centipoise 
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VQ  -   molal  volume  of  the  solute  at  normal   boiling 
point 


Chapter  6 
Results  and  Discussion 

6.1   Equilibrium  studies 

6.1.1  Single  solute  isotherm  studies  for  phenol  and  PCP 

Two  isotherm  experiments  were  conducted,  one  for  phenol 
and  the  other  for  p-chlorophenol .  The  temperature  for  both 
the  experiments  were  maintained  at  25  deg.  C  and  the  pH  was 
maintained  at  7.  Tap  water  was  used  to  make  up  the  solutions. 
This  water  initially  passed  through  a  carbon  column.  three 
feet  in  height  to  and  was  used  to  remove  organics  and  other 
adsorbables  present  in  tap  water. 

The  isotherm  data  was  represented  by  the  three  parameter 
isotherm  equation  shown  below 


q  =  AC/(  1  *  BC   )     for  #  {  1 


This  equation  was  used  for  describing  the  data  for  both 
phenol  and  p-chlorophenol.  The  experimental  and  predicted 
points  are  shown  in  Figures  6-1  and  6-2  for  phenol  and  p- 
chlorophenol,  respectively.  The  plots  are  log-log  plots.  The 
estimation  was  done  by  using  a  non  linear  parameter 
estimation  computer  program  <  Mathews,  1975  ) .  The  parameters 
A,  B,  and  #  are  computed  by  using  the  principal  axis  method 
and  minimization  of  the  sum  of  squares  of  residuals  between 
experimental  and  predicted  values. 

71 


72 


(W9/S310W)    -3N00    aovjyns 


73 


VJ 

Ld 


U'o  z  g  i 

-   O   f  H 


U  3 


'o 


O 


-3  Q 


2 


o 

a   — 
'O 


W9/S310W  '0NO0  HOVdHHS 


74 

From  the  plots  we  see  that  the  curve  for  the  p- 
chlorophenol  isotherm  is  higher  than  that  of  the  phenol 
isotherm.  That  means  that  the  rate  of  adsorption  for  p- 
chlorophenol  is  higher  than  that  of  phenol  adsorption. 

Table  6.1  gives  the  values  of  the  A,  B,  and  /a  for  the 
solutes  phenol  and  p-chlorophenol  at  25  deg.  C  and  pH  7. 


Table  6.1.  Three  parameter  isotherm  constants 

for  phenol  and  PCP  >  i  ,.,  M ,  *  L  *  1 


Solute 


0 


Phenol         15.11      7.546     0.8685 


PCP  30.06      29.92     0.9466 
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6.2  Single  solute  rate  studies  for  single  adsorbent  sizes 

Rate  studies  were  conducted  for  differnt  sizes  of 
carbon  using  a  temperature  of  25  deg .  C  and  pH  of  7.  The 
solutes  used  were  phenol  and  p-chlorophenol .   Tables  6.2   and 

6.3  show  the  operating  conditions  for  batch  studies  for 
phenol  and  p-chlorophenol,  respectively.  The  surface 
dif f usivities  and  batch  film  transfer  coefficient  for  the 
different  sizes  were  estimated  using  a  computer  program 
written  by  Mathews  (  1975  > .  The  computer  program  needed 
an  initial  estimate  for  film  transfer  coefficient  and 
intraparticle  diffusion  coefficient. 

The  film  transfer  coefficient  can  be  calculated  from 
initial  rate  data  according  to  the  formula  <  Mathews  and 
Weber,  1977,  Larson  and  Tien,  1983  >. 

ln(  Cbi  /  Cbo  >  "  "<<  3«M«Kf*t)/(  ap»/op«V)  ) 

Here, 

ap  is  the  radius  of  the  particle 

CDi   is  the  concentration  of  the  solution  at   a   particular 

time 
Cbo  is  the  initial  concentration  of  the  solution 
M  is  the  mass  of  adsorbent 
t  is  the  time 

pp  is  the  density  of  the  particle 
V  is  the  volume  of  the  solution 
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Table  6.2  Operating  conditions  for  batch  studies 
for  phenol  solute 


Sieve        Initial  Weight  of 

size        solution  carbon 

concentration 


(  M/L  )*10+ 


#12-#14 

2.515 

#14-#16 

2.461 

#16-#18 

2.481 

#18-#20 

2.488 

#20-#25 

2.556 

#25-#30 

2.495 

#30-#35 

2.468 

8 

6.5 

3 

8 

8 

6 

6.5 


Table  6.3  Operating  conditions  for  batch  studies 
for  PCP 


Sieve        Initial  Weight  of 

size         solution  carbon 
concentration 

(  M/L  )*10'f-  gms 


#12-#14  2.582 

#18-#20  2.438 

#20-#25  2.561 

#30-#35  2.136 
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The  above  equation  can  be  represented  as  a  linear  graph  and 
holds  good  for  the  first  five  minutes  of  the  experiment. 
However,  this  depends  on  the  solute  and  the  type  of  adsorbent 
used  . 

Figures  6.3  to  6.5  gives  the  experimental  and  predicted 
plots  for  batch  studies  for  phenol  solution  and  size 
fractions  from  #12-#14  to  #30-#35.  Table  6.4  gives  the  sieve 
size  used,  the  film  transfer  coefficients  and  the 
intraparticle  diffusion  coefficients  estimated  from  the  batch 
rate  data  for  phenol . 

From  studies  conducted  by  Mathews  and  Weber  for 
adsorption  in  slurry  reactors  (1984),  they  concluded  that 
mass  transfer  coefficient  was  proportional  to  N0-75,  where  N 
is  the  impeller  speed.  From  our  studies  with  212  micron 
particle  size  we  found  that  the  mass  transfer  coefficient  is 
proportional  to  N0-80. 

Correction  for  kf  is  needed,  especially  in  the  case  of 
adsorption  with  a  mixture  of  adsorbent  sizes.  This  was 
important  in  our  case  of  study  because  batch  studies  were 
conducted  at  two  different  impeller  speeds.  The  impeller 
speeds  used  were  800  rpm  and  700  rpm .  The  equation  for 
conversion  of  mass  transfer  from  one  rpm  to  another  is  given 
by 

0-80 


b 
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Tables    6.5     shows    the    corrected      mass      transfer 
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coefficients  for  the  particles'  in  the  range  of  #20-#25  to 
#30-#35  for  phenol.  The  mass  transfer  coefficients  were 
converted  from  mass  transfer  coefficients  at  impeller  speed 
of  700  rpm  to  mass  transfer  coefficient  at  S00  rpm . 

It  can  be  seen  from  the  experiments  that  the  calculated 
and  estimated  value  of  film  transfer  coefficient  are  very 
close  and  within  5fc  of  each  other.  It  can  also  be  seen  that 
as  the  particle  size  decreased  to  within  a  certain  range, 
the  film  transfer  coefficient  increased  and  then  remained  a 
constant.  Within  the  range  of  #12-#14  to  #18-#20  the  film 
transfer  coefficient  increased  and  within  the  range  of  #20- 
#25  to  #30-#35  the  film  transfer  coefficient  remained  a 
constant  of  almost  around   7.4*10"^  cm/sec. 

It  was  also  seen  that  when  geometric  mean  diameter  was 
used  the  intraparticle  dif f usivities  ranged  from  3.2»10"s 
cm2/sec  to  3.4»10~s  cm2/sec.  but  when  Sauter  mean  diameter 
was  used  the  intraparticle  dif fusivities  remained  a  constant 
at  3.5»10"8  cm2/sec  . 

To  study  the  above  phenomenon  batch  studies  were 
conducted  using  p-chlorophenol  as  the  solute.  Once  agin  it 
was  seen  that  the  intraparticle  diffusion  coefficients 
remained  a  constant  at  3.7«10"8  cm2/sec  when  Sauter  mean 
diameters  were  used  and  ranged  from  2.9»10"s  cm2/sec  to 
3.2«10"°  cm-/sec  when  geometric  mean  diameter  was  used.  The 
film  transfer  coefficients  were  found  to  increase  in  the  same 
range   as  that  in  the  case  of  phenol  and  remained  a   constant 
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Table  6.5  Values  for  film  transfer  coefficient  corrected 
for  impeller  speeds  for  phenol  solute 
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Sieve      Film  transfer  coefficient 
size       for  impeller  rpm  700 
based  on 


Corrected  film 
transfer  coefficient 
for  impeller  rpm  800 
based  on 


Geometric   Sauter 
diameter    diameter 
cm/ sec      cm/ sec 


Geometric     Sauter 
diameter      diameter 
cm/ sec       cm/ sec 


#20-#25  6.26*10~3  7.49*10-3 
#23-#30  6.21*10~3  7.39*10~3 
#30-#35    6.22*10~3   7.42*10-3 


6.96*10-3     8.33*10~3 


6.91*10 


-3 


8.22*10 


6.92*10-3     8.25*10-3 


86 
at  for  the  other  range. 

Figure  6.6  shows  the  plots  of  the  experimental  and 
predicted  values  using  Sauter  mean  as  diameter  using  p- 
chlorophenol  as  the  solute.  Table  6.6  gives  the  values  for  kf 
corrected  for  impeller  speeds  for  p-chlorophenol  adsorption. 
Table  6.7  gives  the  sieve  size  used  and  the  estimated  film 
transfer  coefficients  and  the  intraparticle  diffusion 
coefficients . 
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Table  6.6  Values  for  film  transfer  coefficient  corrected 
for  impeller  speeds  for  PCP  solute. 
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Sieve      Film  transfer  coefficient 
size       for  impeller  rpm  700 
based  on 


Corrected  film 
transfer  coefficient 
for  impeller  rpm  800 
based  on 


Geometric 
diameter 


Sauter 
diameter 


Geometric    Sauter 
diameter     diameter 


cm/ sec 


cm/sec 


cm/ sec 


cm/ sec 


#20-#25 
#30-#35 


9.13*10 


9.20*10 


-3 


10.92*10 


10.93*10 


-3 


10.17*10~3   12.15*10~3 
10.25*10~3   12.16*10~3 
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6.3    Single   solute  rate  studies  for  mixtures   of   adsorbent 
sizes 

Experiments  were  conducted  for  a  mixture  fo  sizes  to 
study  their  effects  on  mass  transfer  and  intraparticle 
diffusion  coefficients  estimated  from  batch  reactor  studies. 
Tables  6.8  and  6.9  gives  the  weight  of  carbon  taken  for  batch 
studies  for  a  mixture  of  particle  sizes  for  phenol  and  p- 
chlorophenol  respectively.  Tables  6.10  and  6.11  give  the 
operating  conditions  for  batch  studies  for  a  mixture  of  sizes 
for  phenol  and  p-chlorophenol  respectively.  Figures  6.7  to 
6.12  show  the  plots  of  the  experimental  and  predicted  curves 
using  Sauter  mean  diameter  for  phenol  adsorption.  The  model 
predictions  were  obtained  using  the  computer  program  for 
single  solute  rote  adsorption  developed  by  Mathews   (1975) . 

When  experiments  were  conducted  for  a  mixture  of 
adsorbent  sizes  some  of  the  individual  size  ranges  fell  in 
the  impeller  rpm  range  of  800  while  the  others  fell  in  the 
range  of  700  rpm.  For  such  cases  the  mass  transfer 
coefficients  for  the  particles  in  the  range  of  700  rpm  were 
corrected  to  mass  transfer  coefficients  for  800  rpm.  The  mass 
transfer  coefficient  for  the  mixture  was  then  calculated  by 
multiplying  the  weight  fraction  of  individual  sizes  present 
in  the  mixture  with  their  respective  mass  transfer 
coefficients  (  corrected  or  uncorrected  depending  on  the 
impeller  speed  used  for  the  experiment  )  and  then  summing 
them   up.   The  intraparticle  diffusion  coefficients  used   for 
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Table  6.9  Weights  taken  (  gms  )  for  mixture  of  sizes  for 
batch  studies  for  PCP  adsorption. 


Mix       Sieve  sizes $> 

no. 

#12-#14   #18-#20   #20-#25  #25-#30 

1  2        2        2         2 

2  0        2        2         4 

3  4        2        2         0 
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Table  6.10  Operating  conditions  for  batch  studies  for 
mixture  of  sizes  for  phenol  adsorption 


Mix  Initial  solution  Weight  of  carbon 

no.  concentration 

(  M  )*103  gms 


1  2.563  8.00 

2  2.522  8.00 

3  2.584  8.00 

4  2.536  8.00 

5  2.386  8.00 

6  2.447  8.00 

7  2.638  8.00 
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Table  6.11  Operating  conditions  for  batch  studies  for 
mixture  of  sizes  for  PCP  adsorption. 


Weight  of  carbon 


Mix 

Initial  solution 

no. 

concentration 

(  M  )*103 

1 

2.383 

2 

2.585 

3 

2.432 

gms 


8.00 
8.00 
8.00 
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geometric  mean  diameters  were  just  the  mean  of  the  individual 
intraparticle  diffusion  coefficients. 

Table  6.12  gives  the  mass  transfer  coefficients  and 
intraparticle  diffusion  coefficients  used  to  predict  the 
kinetic  curves  when  a  mixture  of  particle  sizes  were  used  for 
phenol  adsorption.  From  the  Table  6.12  we  see  that  the 
percent  deviations  for  Mixtures  for  1,  2,  6,  and  7  were  much 
lower  than  percent  deviations  for  Mixtures  3  and  4.  A 
possible  reason  why  Mixtures  1  and  2  are  much  lower  than 
mixtures  3  and  4  could  be  because  the  impeller  speed  used  for 
mixes  1  and  2  was  700  rpm .  At  700  rpm  the  probability  of 
shearing  of  particles  due  to  contact  with  the  impeller  blades 
is  lower.  Moreover.  for  the  Mixtures  3  an  4  the  particles 
were  in  the  range  of  #12-#14  to  #18-#20.  and  in  this  range 
there  was  a  great  deal  of  variation  in  the  mass  transfer 
coefficients  and  the  method  adopted  by  us  for  averaging  these 
mass  transfer  coefficients  might  not  have  been  very  accurate. 

It  is  also  seen  from  Table  6.12  that  the  percent 
deviations  between  experimental  and  predicted  results  were 
much  lower  when  Sauter  mean  of  Sauter  mean  diameters  were 
used  than  when  Sauter  mean  of  geometric  mean  diameters  were 
used  . 

Batch  studies  were  conducted  for  mixture  of  adsorbent 
sizes  using  p-chlorophenol  as  soiute.  In  this  case  too  it  was 
seen  that  Sauter  mean  diameter  of  Sauter  mean  diameters  gave 
better   predictions,    ie.,    the   percent   deviations    from 
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experimental  and  predicted  were  much  lower  when  Sauter  mean 
diameters  of  Sauter  mean  diameters  were  used  than  when  Sauter 
mean  diameters  of  geometric  mean  diameters  were  used.  Figures 
6.13  to  6.15  show  the  experimental  and  predicted  curves 
using  Sauter  mean  diameters  for  p-chlorophenol  adsorption. 
Table  6.13  gives  the  results  of  our  studies  conducted  for  a 
mixture  of  particle  sizes  using  p-chlorophenol  as  the  solute. 
Here  we  see,  again,  that  using  Sauter  mean  of  Sauter  mean 
diameters  gave  better  predictions  than  when  Sauter  mean 
diameters  of  geometric  mean  diameters. 

Hence,  it  was  concluded  that  Sauter  mean  diameters  could 
be  a  better  way  of  representing  particle  diameters  for  batch 
studies  than  geometric  mean  diameter  of  sieve  sizes. 


\&l 


CO 

o 

X 

LjJ 

H 


b  u 

00 
U     3 


&3 

o  -a 

to  <u 

■a  u 

<  a 


J0/0 


105 


q 


o 

. 

m 

rO 

00 

w 

Z> 

CN      ^4 

O 

0) 
•    U 

T 

0   <d 

c  s 

* — " 

a 

■H   -O 

UJ 

S« 

> 

o  <u 

O 

K-t    s 

C\J 

f- 

—   i-. 

25 


&3 

01     <U 

<  a 


5= 

SB 


3/3 


106 


CO 

cr 
o 

X 

UJ 


u 

rH 


■H   T) 

S 


a  c 

O     (0 
M    3 


a  o 
o  ~a 


3 

fa 


O/D 


10? 


nj    O 

•h    a 

T3 

a 

eg   o 
H  -h 

u  CO 

•H  3 

■U  M-i 

M  m      • 

co  -h  a 
a-a  o 

•H 
CD     CU    -U 

oo  cj   a 
(0    co    *4 

a)   j-«   co 

>    3T3 


I 


CO 

•H 

U   *T3     00 
01     4>     C 

> 

eu  ii  *h 

0> 

>!    cj    co 

13 

O-H     3 

CO 

•a 

m 

Jj 

c  at   en 

01 

a 

01    J-i    w 

■■-> 

01 

ci   ari 

a) 

u 

3           3 

■ 

u 

W    T3     CO 

CO 

at   a)   c   aj  -h 


H  T3  w 


•a  ^ 
I 


a)  ,n 

e 

£ 

o 

u    S 

0 

T3 

■a  u 

111 

a  «-> 

CO 

M 

V)  CO 

•O 

a  .a 

co   aj 

u  u 

aj    eg 

3     S 

s 

co  i3  a 

(U     HI    rl 

S     CO     4J 

CO 

co   c: 

•H 

CO 

H   43     CO 

CO 

3 

0)           3 

0)     >» 

O   --* 

4J        •    C 

00  rH 

M 

3    CO 

CO    co 

U    =Q 

CO  -H     fl 

5  a 

•H 

Wfl    O 

■H    Cfl 

S  w 

01 

0) 

U 

■ 

CO 

<U 

CO 

S  »W 

u 

-O 

0    0 

01 

00  c 

3 

03 

CO 

en  -d   o    S    co         6  ^ 


a  c 


108 
6.4   Single  solute  column  studies  for  single  sizes 

The  column  studies  were  the  fixed  bed  type  and  the 
column  was  operated  in  the  downflow  direction.  The  column  was 
operated  with  solutes  being  phenol  and  p-chlorophenol  and 
experiments  were  conducted  to  judge  the  efficiency  of  the 
mathematical  model . 

Table  6.14  gives  the  fixed  bed  operating  conditions  for 
phenol  adsorption.  Figures  6.16  to  6.21  give  the  experimental 
and  predicted  breakthrough  curves  for  phenol  with  single 
sizes.  As  mentioned  earlier,  the  intraparticle  diffusion  rate 
from  batch  studies  were  used,  but  the  external  mass  transfer 
coefficients  used  were  estimated  using  the  method  suggested 
by  Divivedi  and  Upadhyay  <  1979  )  as  discussed  in  chapter  5. 

Table  6,15  gives  the  percent  deviations  between 
experimental  and  predicted  values  along  with  the  diameters  of 
the  particles  and  their  respective  film  transfer  coefficients 
and  intraparticle  diffusion  coefficients. 

It  can  be  seen  from  Table  6.15  that  the  predictions  for 
the  breakthrough  curves  are  closer  to  the  experimental  data 
when  Sauter  mean  diameters  were  used.  It  was  also  seen  that 
there  was  not  much  of  a  pronounced  difference  from  the 
geometric  mean  diameter  and  Sauter  mean  diameter  predictions 
with  respect  to  percent  deviations  from  the  experimental 
points  in  the  case  of  phenol  adsorption. 

To  test  this  the  column  experiments  were  conducted 
again,  only  in  these  cases  using    p-chlorophenol.  Table  6.16 
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Table  6.14  Fixed  bed  operating  conditions  for  phenol  adsorption 


Sieve 
size 

Bed 

diameter 

Bed 
height 

Flow 
rate 

Hydraulic 
loading  rate 

Weight  of 
carbon 

cm 

cm 

ml/min 

m/min 

gms 

)  #12-#14 

5.08 

30.9 

500 

0.247 

250 

■  #14-#16 

5.08 

30.3 

500 

0.247 

250 

#16-#18 

5.08 

29.9 

500 

0.247 

250 

#18-#20 

5.08 

29.1 

500 

0.247 

250 

■  #20-#25 

3.1 

23.2 

125 

0.166 

75 

#25-#30 

3.1 

22.9 

125 

0.166 

75 

#30-#35 

3.1 

22.8 

125 

0.166 

75 
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Table  6.16  Fixed  bed  operating  conditions  for  PCP  adsorption 


Sieve     Bed         Bed       Flow      Hydraulic      Weight  of 
size      diameter     height    rate      loading  rate   carbon 


ml/min    m/min  gins 


#12-#14 

5.08 

30.9 

500 

0.247 

250 

#18-#20 

5.08 

29.1 

500 

0.247 

250 

#20-#25 

3.1 

23.2 

125 

0.166 

75 

#30-#35 

3.1 

22.8 

125 

0.166 

75 
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shows  the  fixed  bed  operating  conditions  for  p-chlorophenol 
adsorption.  Figures  6.22  and  6.23  shows  the  experimental  and 
predicted  breakthrough  curves  for  p-chlorophenol.  Table  6.17 
shows  the  percent  deviations  for  the  predictions  from 
experimental  points  when  geometric  mean  diameters  and  Sauter 
mean  diameters  were  used.  It  was  seen  once  again  that  Sauter 
mean  diameters  gave  better  predictions  than  when  geometric 
mean  diameters  were  used,  but,  in  this  case,  the  predictions 
were  more  pronounced  as  opposed  to  that  for  phenol  fixed  bed 
predictions . 
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S.5   Single  solute  column  studies  for  mixture  of  carbon  sizes 

Column  studies  were  conducted  for  mixture  of  carbon 
sizes.  Tables  6.18  and  6.19  give  the  weights  of  carbon  sizes 
taken  for  column  studies  for  phenol  and  p-chlorophenol 
experiments,  respectively.  Table  6.20  gives  fixed  bed 
operating  conditions  for  phenol  adsorption  for  mixture  of 
sizes.  Figures  6.24  to  6.26  give  the  predictions  and 
experimental  breakthrough  curves  for  the  mixture  of  sizes 
with  phenol  as  the  solute.  Table  6.21  gives  the  percent 
deviation  for  the  predictions  from  the  experimental  and  the 
diameters  and  their  respective  film  transfer  coefficient  and 
intraparticle  diffusion  coefficient.  In  these  cases  the  mean 
diameters  were  determined  first.  Then,  using  the  equation  by 
Upadyay  and  Dwiwedi  (1977).  the  film  transfer  coefficients 
were  determined.  Surface  diffusion  coefficients  were 
determined  by  taking  the  mean  of  the  surface  diffusion 
coefficients  of  the  individual  particles  which  were 
established  from  batch  studies  with  single  sized  particles. 

From  the  experiments  it  was  seen  that  Sauter  mean  of  the 
Sauter  mean  diameters  gave  better  predictions  than  when 
Sauter  mean  of  geometric  mean  diameters  were  used.  This  same 
occurrence  was  noticed  in  the  case  of  p-chlorophenol 
adsorption  for  a  mixture  of  sizes.  Table  6.22  gives  the  fixed 
bed  operating  conditions  for  p-chlorophenol  adsorption. 
Figures   6.27   and  6.23  give   the  breakthrough  curves   for   a 
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Table  6.19  Weights  taken  ("gins  )  for  mixture  of  sizes  for 
fixed  bed  studies  for  PCP  adsorption. 


Mix     Sieve  sizes > 

no. 

#12-#14   #18-#20   #20-#25   «0-#35 


1 

60 

60 

60 

70 

2 

80 

80 

0 

90 

3 

0 

0 

37.5 

37.5 

125 


Table  6.20  Fixed  bed  operating  conditions  for  phenol  adsorption 
for  a  mixture  of  sizes 


Mix 

no. 

Bed 
diameter 

Bed 

height 

Flow 
rate 

Hydraulic 
loading  rate 

Weight  of 
carbon 

cm 

cm 

ml/mln 

m/min 

gms 

1 

5.08 

30 

500 

0.247 

250 

2 

5.08 

28.7 

500 

0.247 

250 

3 

3.1 

22.9 

125 

0.166 

75 
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5.08 

29.6 

500 

0.247 

250 
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Table  6.22  Fixed  bed  operating  conditions  for  PCP  adsorption 
for  mixture  of  sizes. 


Mix 
no. 

Bed 
diameter 

Bed 

height 

Flow 
rate 

Hydraulic 
loading  rate 

Weight  of 
carbon 

cm 

cm 

ml/min 

m/min 

gms 

1 

5.08 

26.5 

500 

0.247 

250 

2 

5.08 

27.6 

500 

0.247 

250 

3 

3.1 

23.0 

125 

0.166 

75 
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mixture  of  sizes  for  p-chlorphenol  adsorption.  Table  6.23 
gives  the  percent  deviation  of  predicted  from  experiment ai 
for  mixture  of  sizes  p-chlorophenoi  fixed  bed  operations. 
From  the  experiments  it  was  noticed  that  Sauter  mean 
diameters  gave  better  predictions  than  when  geometric  mean 
diameters  were  used. 

The  computer  program  was  also  run  taking  into 
consideration  a  layered  bed.  The  individual  film  transfer 
coefficients  and  the  individual  surface  diffusion 
coefficients  replaced  the  averaged  film  transfer  coefficients 
and  surface  diffusion  coefficients  and  individual  particle 
diameters  replaced  averaged  particle  diametrs.  Using  this 
approach  it  is  seen  in  Tables  6.21  and  6.23  that  the  percent 
deviations  between  experimental  and  predicted  data  are  much 
closer  than  either  Sauter  mean  of  Sauter  mean  diameters  or 
Sauter   mean   of  geometric  mean  diameters. 

From  fixed  bed  studies  we  concluded  that  the  Sauter 
mean  diameters  was  a  better  way  of  representing  particle  size 
distribution   for  column  studies. 
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6.6      Multicomponent    adsorption 

6.6.1  Experimental  methods 

Experimental  procedures  for  multicomponent  adsorption 
were  similar  to  that  of  single  solute  adsorption,  i.e.,  all 
the  preliminary  precautions  taken  for  single  solute 
adsorption  were  applicable  for  multicomponent  adsorption.  In 
all  studies  (  equilibrium,  batch  and  fixed-bed  studies  )  the 
temperature  was  maintained  at  25  deg .  C  and  pH  was  maintained 
at  7. 

6.6.2  Equilibrium  studies 

In  multicomponent  equilibrium  studies,  experimental  data 
for  phenol  and  p-chlorophenol  was  collected.  To  describe  the 
data  two  multicomponent  models  were  selected.  They  were 

1)  Mathew's  model  (1975)  which  is  given  as 


%»A&hd/0  +  Bl&lrlt) 


2)  Prausnitz's  model  <  IAS  model  ) . 
The. isotherm  constants  A,   B,   and  #   used  for  the  models  were 
obtained  from  single  solute  isotherm  studies. 

Recently,  it  has  become  necessary  to  make  use  of 
theoretical  models  that  allow  the  predictions  of 
multicomponent  adsorption  with  less  amount  of  experimental 
data.  Tables  6.24  and  6.25  show  the  experimental  and 
predicted  values  for  the  multicomponent  studies  using  the 
models  stated  above.   From  the  tables  we  see  that  both  modeis 
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gave  good  predictions  for  the  experimental  data.  It  was  seen 
that  Wathew's  model  was  under  predicting  for  high 
concentrations.  The  values  of  Tft  for  ilathew's  model  was 
0.4679  and  0.1945.  The  percent  deviations  for  phenol  solute 
and  p-chlorophenol  solute  was  estimated  using  the  equation  as 
follows : 

«         V       W  N       ' 


Using  the  above  equation  it  was  seen  for  the  Mathew'a  model 
the  percent  deviation  between  experimental  and  predicted 
values  for  surface  adsorption  was  24*  whereas  for  the  IAS 
model  the  percent  deviation  was  14*. 
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6.6.3  Batch  studies 

Multicomponent  adsorption  for  batch  studies  were 
conducted.  Batch  experiments  were  conducted  for  two  single 
adsorbent  sizes  and  one  mixture  of  adsorbent  sizes.  The  aim 
in  these  experiments  was  to  see  how  the  multicomponent 
adsorption  rate  model  would  predict  batch  data.  The  IAS 
equlibrum  model  was  incorporated  into  the  multicomponent  rate 
model . 

Figures  6.29  and  6.30  gives  the  experimental  and 
predicted  values  for  the  multicomponent  adsorption  rate  model 
written  by  Mathews  (1983).  The  program  needed  inputs  for 
film  transfer  coefficients  and  surface  diffusion 
coefficients.  It  also  needed  the  constants  of  the  three 
parameter  isotherm  equation.  All  these  inputs  were  determined 
from  single  solute  rate  and  equilibrium  studies.  The 
diameters  used  in  these  predictions  were  the  Sauter  mean 
diameter  and  the  geometric  mean  diameter.  Table  6.26  gives 
the  operating  conditions  for  batch  studies  for  multicomponent 
adsorption.  Table  6.27  gives  the  percent  deviations  from 
predicted  and  experimental  values  for  multicomponent  rate 
studies.  From  this  table  we  see  that  Sauter  mean  diameters 
gave  better  predictions  than  when  geometric  mean  diameters 
were  used.  It  was  also  seen  that  the  incorporation  of  the  IAS 
model  for  multicomponent  rate  studies  gave  predictions  which 
were  very  close  to  the  experimental  data. 

Batch    studies   for   multicomponent    adsorption    were 
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Table  6.27  Percent  deviations  between  experimental 
and  predicted  results  for  bisolute  batch 
studies  using  diameters  geometric  and 
Sauter. 
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Sieve 

%  deviation 

%  deviation 

size 

for  geometric 

for  Sauter 

mean 

mean 

#12-#14 

8.23 

4.1 

#18-#20 

7.31 

0.99 

Mix  4 

1.27 

4.29 

m 

conducted  for  a  mixture  of  carbon  of  sizes.  The  film  transfer 
coefficients  and  the  surface  diffusion  coefficient  were 
determined  the  same  way  as  done  in  the  case  of  single  solute 
batch  studies  for  a  mixture  of  adsorbent  sizes.  Figure  6.31 
gives  the  experimental  and  predicted  data  for  the  mixture  of 
solutes  using  a  mixture  of  adsorbents.  Table  6.26  shows  the 
percent  deviation  between  experimental  and  predicted  values 
for  Sauter  mean  diameter  and  geometric  mean  diameter. 

6.6.4  Fixed-bed  studies 

A  total  of  four  column  studies  were  conducted,  three  for 
single  adsorbent  sizes  and  one  for  a  mixture  of  adsorbent 
sizes.  No  predictions  were  conducted  and  the  Figures  6.32  to 
6.34  show  the  experimental  values  for  the  column  runs  for 
single  adsorbent  sizes. 

From  these  plots  we  see  that  as  the  particle  3ize 
decreased  in  diameter  the  peak  values  attained  for  phenol 
concentration  in  the  effluent  increased.  It  was  also  seen 
that  C/C0  values  in  all  of  these  runs  exceeded  1.0.  The 
possible  reason  for  this  could  be  due  to  the  fact  that  p- 
chlorphenol  is  adsobed  faster  at  the  top  of  the  column 
initially  and  phenol  at  the  bottom  of  the  column.  As  the  top 
gets  saturated  p-chlorophenol  moves  downward  displacing  the 
adsorbed  phenol  back  into  the  solution  and  this  raised  the 
concentration  of  phenol  in  the  water. 

Figure   6.35   shows   the  experimental   data   for   column 
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studies  for  a  mixture  of  carbon  sizes.  A  total  of  250  gms  of 
carbon  was  taken.  From  the  plot  we  see  that  the  peak  value 
attained  for  phenol  solute  in  the  effluent  was  greater  than 
the  peak  value  attained  for  single  adsorbent  size  #12-#14, 
but  was  lower  than  that  for  #18-#20  and  *>30-#35.  For  sieve 
#03-#35  the  weight  of  carbon  taken  was  75  gm . 

The   initial  concentrations  for  the  influents  for   sieve 
sizes  #12-#14,  #18-#20.  and  #30-#35  was  1 . 847«10-3»  J  (  ' 
1.643»10"3.  and  1.900»10~3  for  phenol  concentrations.  POP 
concentrations  were  1.831*10~3,  1.751»10~3.  and 
1.817»10~3.    For   mixture   of   sizes   the   initial   influent 
concentrations  was  1.832»10~3  for  phenol,  and  for  PCP  it  was 
1.862«10"3. 
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Chapter  7 
Summary  and  Conclusions 


Equilibrium  studies  were  conducted  for  the  adsorption  of 
two  solutes  1)  phenol  and  2)  parachlorophenol  on  activated 
carbon  from  water.  The  tests  were  conducted  for  both  solutes 
at  25  deg.  C.  and  the  pH  of  the  tapwater  was  adjusted  to  7. 
(  The  water  initially  passed  through  a  three  foot  carbon 
column  to  remove  organic  and  other  adsorbables  from  it.  Both 
solutes  took  three  days  to  attain  equilibrium  ) .  The  data 
were  fitted  by  the  three  parameter  isotherm  equation  which 
was  established  by  Mathews  and  Weber  <1977).  It  was  seen  that 
the  three  parameter  isotherm  fitted  the  experimental  data 
well.  In  the  case  of  p-chlorophenol  the  adsorption  isotherm 
curve  was  higher  than  for  phenol  isotherm  indicating  that  p- 
chlorophenol  is  more  strongly  adsorbed  than  phenol. 

Batch  studies  were  conducted  for  phenol  and  p- 
chlorophenol  for  single  adsorbent  sizes  as  well  as  for  a 
mixture  of  carbon  sizes.  Two  different  methods  of  averaging 
the  particle  sizes  were  tried.  They  are  the  geometric  mean 
diameters  and  Sauter  mean  diameters.  For  all  the  experiments 
conducted  during  batch  studies  the  pH  of  the  water  was 
maintained  at  7  and  the  temperature  regulated  to  25  deg.  C. 

For  batch  studies  involving  phenol  and  p-chlorophenol 
adsorption  and   single  sizes,   the  film  transfer   coefficient 
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increased  as  particle  aiza  decreased  to  within  the  size  range 
of  1544  microns  to  917  microns.  It  then  remained  constant 
within  the  size  range  of  772  microns  to  543  microns. 

It  was  also  seen  that  when  Sauter  mean  diameter  based  on 
area  was  used  the  surface  diffusion  coefficient  remained  a 
constant  for  both  phenol  and  p-chlorophenol  at  3.5*10~6  and 
3.7*10~s  cm2/sec.  respectively.  Whereas,  when  geometric  mean 
diameters  were  used  the  surface  diffusion  coefficient  varied 
from  3.2*10~a  to  3.4»10"a  for  phenol  solutes  and  varied  from 
2.9»10"s  to  3.2»10"8  for  p-chlorophenol  solute. 

In  the  case  of  batch  studies  for  a  mixture  of  sizes,  it 
was  seen  that  better  predictions  for  the  data  were  obtained 
when  the  Sauter  mean  diameters  were  used  than  when  geometric 
mean  diameters  were  used  to  average  the  particle  size 
distribution.  For  both  phenol  and  p-chlorophenol  batch 
studies  for  a  mixture  of  adsorbent  sizes,  it  was  seen  that 
the  predicted  and  experimental  data  were  within  5fc  of  each 
other . 

Column  studies  were  conducted  for  single  adsorbent  sizes 
as  well  as  for  a  mixture  of  adsorbent  sizes.  Columns  of  two 
different  diameters  were  used  depending  on  the  size  of  the 
adsorbent.  The  temperature  and  pH  were  maintained  at  25  deg. 
C  and  7,  respectively.  The  program  used  to  predict  the 
breakthrough  profiles  was  the  homogenous  solid  phase 
diffusion  model  which  was  modified  by  Mathews  <1984')  so 
that  equilibrium  could  be  represented  by  the  three  parameter 
isotherm . 
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For  column  studies  it  was  seen  that  when  Sauter  mean 
diameter  was  used,  the  model  predictions  were  much  closer  to 
the  experimental  data  than  when  geometric  mean  diameter  was 
used.  This  was  noticed  in  the  case  of  single  adsorbent  sizes 
as  well  as  for  a  mixture  of  adsorbent  sizes.  In  both  cases, 
ie.,  for  single  adsorbent  sizes  as  well  as  for  a  mixture  of 
adsorbent  sizes  the  percent  deviations  from  experimental  and 
predicted  data  were  within  5fc  of  each  other. 

Bulticomponent  equilibrium  studies  were  conducted  in  a 
similar  fashion  as  that  of  the  single  solute  equilibrium 
studies.  Two  models  were  used  to  fit  the  isotherm  data.  They 
were  the  extension  of  the  three  parameter  isotherm  equation 
established  by  Bathew's  <  1977  )  and  the  ideally  adsorbed 
solution  model  (  IAS  model  ) .  However,  the  percent  deviations 
between  experimental  and  predicted  data  were  lower  when  the 
IAS  model  was  used  than  when  Bathew's  model  was  used.  When 
using  Bathew's  model  the  percent  deviation  was  24*  whereas 
when  using  the  IAS  model  the  percent  deviation  was  14*. 

Batch  studies  were  conducted  for  single  as  well  as  for  a 
mixture  of  sizes.  Predictions  were  obtained  by  making  use  of 
the  multicomponent  adsorption  rate  model  program  written  by 
Bathew's  <  1983  ).  It  was  seen  that  the  use  of  Sauter  mean 
diameters,  based  on  area,  gave  better  predictions  than  when 
geometric  mean  diameters  were  used.  This  was  also  noticed  in 
the  case  of  batch  studies  for  a  mixture  of  adsorbent  sizes. 

In  the  case  of  fixed-bed  studies  with  mul ticomponent3  it 
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was  seen  that  as  the  particle  size  decreased  the  peak  value 
attained  for  phenol  solute  in  the  effluent  kept  increasing 
with  time.  It  was  also  noticed  that  C/C0  values  for  phenol 
adsorption  exceeded  1.0.  The  possible  reason  for  this  could 
be  due  to  the  fact  that  p-chlorophenol  is  adsorbed  faster  at 
the  top  of  the  column  initially  and  phenol  at  the  bottom  of 
the  column.  As  the  top  became  saturated,  p-chlorophenol  moved 
downward  displacing  the  adsorbed  phenol  back  into  the 
solution . 

For  fixed  bed  studies  for  multicomponent  adsorption 
using  mixture  of  adsorbent  sizes,  it  was  seen  that  the  peak 
value  attained  for  phenol  solute  in  the  effluent  was  higher 
than  the  peak  value  attained  for  multicomponent  adsorption 
for  single  adsorbent  size  #12-#14,  but  was  lower  than  that 
for  #18-#20  or  #30-#35. 


APPENDIX 


155 


156 


00 
0) 

K 

1 

LU 

<M 

1- 
2 

=** 

o 

< 

£ 

a 

Q 

N- 

X 

Q 
LU 

1— 

< 

h. 

H 
X 
<£ 

LU 

3 

lu 
> 


2 
o 


< 

X 

I- 

<E 
LU 


1    1    I  -I- 

O  0)  CD    Is- 

1 

10 

1 
in 

i 

T 

— 1 

-  o  o  o 

o 

O 

o 

O 

(iuo)    M3X3^Via 


(  Bra  )  aaisiwia 


en 

__  en 


1  "S 


co 

'  0) 


o 

0) 

o 

o 

CM 


- 1- 


O01CD 

-  6  6 


~r- 


in 
O 


— P" 


ro 

o 


(  <no  )  aaiaHvia 


(  ra3  ) 


OfflO)    K 

10 

in 

T 

ro 

-  o  o  o 

o 

O 

o 

O 

aaiaiwia 

1<0 


(  ma  )  naiaHVia 


1:1 


o 

•w   en 

fO 

3     QJ 

tH 

<u    a 

>  -H 

<H   -u 

4-1      J-t 

CO    (3 

(   "13    )    H3I3HVICI 


152 


n  in 

U  H 


(  M3  )  aaaawia 


153 


_o 


— r- 
O 


O 


— i 1 — 

o      a 
—      6 


o 

H 

J-l 

1^ 

aj 

h 

H 

C) 

•H 

■U 

M 

n> 

V 

H 

a. 

o 

•H 

a 

lO 

T3 

o 

ID 

d 


d 


(  Q3  )  raiawvia 


3  a 

a 

o  vD 
•H  rH 
U  ^ 
U     I 


S   u 


l&i 


U     U 


OS    M 

0)     tU 

o 

iH    U 

h- 

m  a 

OJ    -H 

W     l-i 

OJ     V 

S   a 

a 

O 

•H     C 
13     O 

ID 

jd  tj 

4-1     QJ 

•H      CO 

3    rt 

ja 

o 

O   CO 

fO 

■H    H 

(  °D  )  Haiawvic 


165 


T3  -H 

CU  JJ 

U  u 

01  M 

u  a 

CO 
01 

a  j= 

cc  JJ 

4J  U-l 


-U     U 

6   o. 


o  o 

■H  cn 

O  I 

Bj  SO 

M  rH 

U  CU 


(  ma   )   H3I3KVia 


16< 


1  s. 


-O 

■H    CS 

4J    =S; 

ro 

O       1 

M  cN 

00  -H 

at 


,00 
0) 


o 


O 
<0 


i<? 


<D    A 

4-1 

OJ     01 

■  -c 

o 

■H 

CM 

T3   <W 

O 

J3 

—f— 


'    i    i    T" 

Offl  CO    N 

1 

1 

in 

— r- 

rn 

-  o  o  o 

o 

o 

o 

O 

(  mo  )  aaiawvia 


Si 

a 

— 1 

'} 

t 

::j 

<n 

V 

T, 

> 

£ 

-H 

o 

fl 

rO 

H 

-l 

u"l 

3 

(N 

U 

* 

168 


05 
0) 


CO 


o 

0) 


o 


4J     U 

I1 


0 

u"l 

•d 

m 

u 

o 

u 
en 

i 

CVJ 

M 

pi 

I- 1 1 1- 

Q01CO    s 

— I— 
to 

1 

ID 

<r 

ro 

—  o  o  o 

o 

o 

o 

O 

(  mo  )  raiawvici 


References 


Allen  T.  (1973).  Particle  Size  Measurment  2nd  edition. 

Butler,  J.  A.  V.  and  C.  Ockrent  (1930).  Studies  in 
Electrolocapillarity .  3  J.  Phys .  Chem.  34,2841. 

Cooney  D.  0.  and  F.  P.  Strusi  (1972).  Analytical  Description 
of  Fixed  Bed  Sorption  of  Two  Langmuir  Solutes  Under 
Nonequilibrium  Conditions.  Ind.  Eng.  Chem.  Fund.,  11,123. 

.Crittenden  J.  C.  (1976).  Mathematical  Modelling  of  Fixed  Bed 
Adsorber  Dynamics  -  Single  Component  and  Multiconponent . 
PH.d.  Thesis,  Univ.  of  Michigan,  Ann  Arbor,  MI. 

J/Crittenden  J.  C.  and  W.  J.  Weber.  Jr.  (1978).  Mathematical 
Modelling  of  Fixed  Bed  Adsorption  Systems:  Single  Component 
and  Multicomponent .  71st  Ann.  Meeting,  flier .  Inst.  Chem. 
Eng.,  Miami  Beach,  Florida. 

Dietz  F.  ,  and  J.  Trand.  Vom  Wasser  51:233-257  (1978). 

-Digiano  F.  A.  and  W.  J.  Weber,  Jr.  (1973).  Sorption  Kinetics 
in  Infinite  Batch  Experiments.  J.  Water  Poll.  Control  Fed., 
45,  713. 

Dwiwedi  P.  N.  and  Upadyay  S.  N.  (1977).  Ind.  Eng.  Chem.. 
Process  Des.  Dev . ,  Vol.  IS,  No. 2,  1977. 

•'Famularo  J.,  J.  A.  Mueller  and  A.  S.  Pannus  (1980). 
Predictions  of  Carbon  Column  Performance  from  Pure-Solute 
Data.  J.  Water  Poll.  Control  Fed.  52,2019. 

*  f  oust  A.  5.,  Wenzel  L.  A.,  Clump  C.  W.,  Maus  L.,  Andersen  L. 
B.  Principles  of  Unit  Operations,  pg  525.  Appendix  B. 

f  Fritz  W..  W.  Merk,  E.  U.  Schlunder  and  H.  Sontheimer  (1980). 
Competitive  Adsorptions  of  Dissolved  Organics  on  Activated 
Carbon.  In:  Activated  Carbon  Adsorption  Of  Organics  From  The 
Aqueous  Phase.;  Vol.  1.  (  M.  J.  McGuire  and  I.  P.  Suffet, 
Eds.  ).  Ann  Arbor  Sci . ,  Ann  Arbor,  MI. 

Gleuckauf,  E.  and  J.  I.  Coates  (1947).  Theory  of 
Chromatography,  Part  4.  J.  Chem.  Soc . ,  1315. 

Hall  K.  R.,  L.  C.  Eagleton,  A.  Acrivos  and  T,  Vermeulen 
(1966).   Pore   and   Solid  Diffusion  in  Fixed  Beds   Adsorption 

159 


170 


Under  Constant  Pattern  Conditions.  Ind.  Eng.  Chem.  Fund. 
5,212. 

/Harold  Hey wood ,   Symposium  on  Particle  Size  Analysis.   pg  14. 
(1947),  Institution  o£  Chemical  Engineers. 

Heister  N.  K.  and  T.  Vermeulen  (1952) .  Saturation  Performance 
of  Ion-Exchange  and  Adsorption  Columns.  Chem.  Eng.  Prog., 
48,505. 

Holzel  G..  F.  Fuchs  and  G.  Baldauf  (1979).  Ubersuchunger  Zur 
Optimierung  der  Akti vkohleanwendung  bei  der 
Trikwasserauf bereitung  am  Rhein  unter  besonderer 
Berucksichtigung  der  Regeration  nach  Thermischen  Verfahven. 
In:  Optimeirung  der  Akti vkohleanwendung  bei  der 
Trikwasserauf bereitung.   Universitat   Karlsruhe   Deutschland. 

•\Hseih  J.  s.  c.,  R.  M.  Turian  and  C.  Tien  (1977).  Multi- 
Component  Liquid  Phase  Adsorption  In  Fixed  Beds.  Amer .  Inst. 
Chem.  Eng.  J.,  23,263. 

Jain  J.  S.,  V.  L.  Snoeyink  (1973).  Adsorption  From  Bisolute 
Systems  on  Activated  Carbon.  J.  Water  Poll.  Control  Fed. 

Jossens  L.,  J.  M.  Prausnitz,  W.  Fritz,  E.  U.  Schlunder  and  A. 
L.  Myers  (1978).  Thermodynamics  of  Multisolute  Adsorption 
From  Dilute  Aqueous  Solutions.,  Chem.  Eng.  Sci . ,  33,1097. 

Keinath  T.  M.  and  W.  J.  Weber,  Jr.  (1968).  A  predective  Model 
for  The  Design  of  Fluid  Bed  Adsorbers.  J.  Water  Poll.  Control 
Fed.  40,741 

Keinath  T.  M.  (1977).  Design  and  Operation  of  Activated 
Carbon  Adsorbers  Used  for  Industrial  Wastewater 
Decontamination.  Amer.  Inst.  Chem.  Eng.  Symp.  Ser .  73,1. 

Kyte  W.  S.  (1973).  Nonlinear  Adsorotion  in  Fixed  Beds.  Chem. 
Eng.  Sci.  28,153 

Langmuir  J.  (1918).  The  Adsorption  of  Gases  on  Planes  of 
Glas3,  Mica,  and  Platinum.  J.  Amer.  Chem.  Soc.  40,1361. 

\  Larson  A.  C.  and  C.  Tien  (1983).   Multicomponent   Liquid 

Phase  Adsorption  in  Batch  .  Part  2  Experiment  on  Carbon 
Adsorption  From  Solutions  of  Phenol,  0-Cresol,  and  2,4- 
dichlorophenol .  Chem.  Eng.  Commun .  27,  359-361 

Lee  M.  C,  J.  C.  Crttenden ,  V.  L.  Snoeyink  and  W.  E.  Thacker 
(1980).  Mathematical  Modelling  Of  Humxc  Substances  Removal 
With  Activated  Carbon  Beds.  National  Conference  on 
Environmental  Engineering.  ASCE . ,  NY. 

<j,  Liapis  A.   I.   and  D.   W.  T>  Rippen  (1978).  The  Simulation  of 


171 

Binary  Adsorption  in  Activated  Carbon  Columns  U3ing  Estimates 
of  Diffu3ional  Resistances  Within  the  Carbon  Particles 
Derived  from  Batch  Studies.  Chen.  Eng.  Sci .  35,593 

j  Mathews  A.  P.  (1975).  Mathematical  Modelling  of 
Multicomponent  Adsorption  in  Batch  Reactors.  PH.D.  Thesis. 
Univ.  of  Michigan,  Ann  Arbor,  MI. 

Mathews  A.  P.  and  W.J.  Weber,  Jr.  (1977).  Effects  of  External 
Mass  Transfer  and  Intraparticle  Diffusion  on  Adsorption  Rates 
in  Slurry  Reactors.  Americ.  Inst.  Chem .  Eng.  Symp.  Ser . , 
73,91. 

Mathews  A.  P.  and  S.  R.  Kulkarni  (1983).  Unpublished  work. 
Department  of  Civil  Engineering,  Kansas  State  University. 

Mathews  A.  P.  and  C.  A.  Su .  (1983).  Prediction  of  Competitive 
Adsorption  Kinetics  for  Two  Priority  Pollutants. 
Environmental  Progress  2(4),  257. 

Mathews  A.  P.  (1983).  Adsorption  In  An  Agitated  Slurry  of 
Polydisperse  Particles.  AIChE  Symposium  Series  79,  18. 

Mathews  A.  P.  (1984).  Unpublished  work.  Department  of  Civil 
Engineering,  Kansas  State  University. 

[  Mathews  A.  P.  and  W.  J.  Weber,  Jr  (1984).  Modelling  and 
Parameter  Evaluation  of  Adsorption  in  Slurry  Reactors.  Chem. 
Eng.  Commun.  25,  157. 

Radke  C.  J.  and  J.  M.  Prausnitz  (1972).  Adsorption  of  Organic 
Solutes  from  Dilute  Aqueous  Solution  on  Activated  Carbon. 
Ind.  Eng.  Chem.  Fund.,  11,445. 

Radke  C.  J.  and  J.  M.  Prausnitz  (1972).  Thermodynamics  of 
Multisolute  Adsorption  from  Dilute  Liquid  Solutions.  Amer. 
Inst.  Chem.  Eng.  J.,  18,761. 

.»  Schay  G.  J.,  F.  P.  Fejes  and  J.  Szathmary  (1957).  Studies  on 
the  Adsorption  of  Gas  Mixture,  1.  Statistical  Theory  of 
Physical  Adsorption  of  Langmuir  type  in  Multicomponent 
Systems.  Acta  Chem.  Sci.  Hungary.,  12,299. 

Snoeyink,  V.L.  and  W.J.  Weber,  Jr.  (1968).  Reaction  of  the 
hydrated  proton  with  Active  Carbon.  In:  Adsorption  from 
Aqueous  Solution,  Advances  in  Chemistry  Series,  No.  79,112. 

Suzuki  M.  and  Kawazoe  (1974).  Batch  Measurment  of  Adsorption 
Rate  in  an  Agitated  Tank  -  Pore  Diffusion  Kinetics  with 
Irreversible  Isotherm.  J.  Chem.  Eng.  Japan,  7,346. 

Sweeny  M.  W.,  W.  A.  Melville,  B.  Trgovicich  and  C.  p.  i_. 
Grady,   Jr.  (1982).  Adsorption  Isotherm  Parameter  Estimation. 


172 
J.  Env.  Eng.  Div.,  ASCE,  108,913. 

Thacker  W.  E.  and  V.  L.  5noeyink  (1978).  On  Evaluating 
Granular  Activated  Carbon  Adsorption.  Notes  and  Comments,  J. 
Amer.  Water  Works  Assoc,  70,45. 

Thomas  H.  C.  (1944).  Hetrogenous  Ion  Exchange  in  a  Flow 
System.  J.  Amer.  Chem .  Soc . ,  66,1664. 

Van  hall  C.  E.,  Ed,  Measurments  of  Organic  Pollutants  in 
Water  and  Wastewater  Treatment,  ASTM  Spec.  Pub  686 
(Philadelphia:  American  Society  for  Testing  and  Materials, 
1979  ) . 

Van  Vliet,  B.  M.  and  W.  J.  Weber,  Jr  (1979).  Comparitive 
Performance  of  Synthetic  Adsorbents  and  Activated  Carbon  for 
Specific  Compound  Removal  from  Wastewaters.  52nd  Ann.  Conf . , 
Water  Poll.  Control  Fed.,  Houston,  TX . 

Van  Vliet,  B.  21.  and  W.  J.  Weber,  Jr.  (1980).,  Modelling  and 
Prediction  of  Specific  Compound  Adsorption  by  Activated 
Carbon.  Water  Research,  14(12), 1719. 

Vassiliou  B.  and  J.  5.  Dranoff  (1962).  The  Kinetics  of  Ion 
Exclusion.  Amer.  Inst.  Chem.  Eng.  J.,  8,248. 

Vermeulen  T.  (1953).  Theory  for  Irreversible  and  Constant 
Pattern  Solid  Diffusion.  Ind.  Eng.  Chem.  Fund.,  45. 

Weber  W.  J.,  Jr.  Physicochemical  Process  For  Water  Quality 
Control.  (1972),  Wiley-Interscience . 

Weber  W.  J..  Jr.  and  R.  R.  Rumner  (1965).  Intraparticle 
Transport  of  Sulphonated  Alkyl  Benzenes  in  a  Porous  Solid: 
Diffusion  with  Nonlinear  Adsorption.  Water  Resources  Res., 
1,361. 

Westermark  M.  (1975).  Kinetics  of  Actvated  Carbon  Adsorption. 
J.  Water  Poll.  Control  Fed.,  47,704. 

Zwiebel  I.,  R.  L.  Gariepy  and  J.  J.  Schnitzer  (1972).  Fixed 
Bed  Desorption  Behavior  of  Gases  with  Non-Linear  Equilibria  : 
Part  1.  Dilute,  One  Component  Isothermal  Systems.  Amer.  Inst. 
Chem.  Eng.  J.,  18,1139. 


EFFECT  OF  PARTICLE  SIZE  DISTRIBUTION 
ON  ACTIVATED  CARBON  ADSORPTION 


by 
THOPPIL  JOJO  KUNJUPALU 


B.E.,  Collage  of  Engineering  Guindy,  Madras  University 
India,  1963 


AN  ABSTRACT  OF  A  WASTER'S  THESIS 


submitted  in  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Civil  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 


1986 


Abstract 

Batch  and  fixed  bed  exper iments  were  conducted  to  study 
the  effect  of  particle  size  on  carbon  adsorption.  Three 
different  types  of  diameters  were  used.  They  were  geometric 
mean  diameter  of  the  sieve  sizes,  length  mean  diameter.  and 
Sauter  mean  diameter.  All  these  diameters  were  estimated 
based  on  area  of  the  particles  using  Quantimet-720  image 
analysis.  A  mathematical  model  was  used  to  compare  the 
results  of  batch  and  fixed  bed  experiments.  The  model 
equations  were  formed  considering  film  transfer  and  surface 
diffusion  controlled  adsorption  rates.  For  fixed  bed  studies 
the  model  equations  were  solved  by  using  orthogonal 
collocation  technique. 

Phenol  and  parachlorophenol  were  the  two  solutes  used  in 
the  experimental  studies.  All  experiments  were  conducted 
using  tapwater  at  25  deg .  C  and  pH  7.  Single  solute  isotherm 
data  were  fitted  by  the  three  parameter  isotherm  equation 
developed  by  Mathews  and  Weber  (1977).  Bisolute  isotherm 
studies  were  conducted  using  Mathews  model  (1975)  and  the  IAS 
model . 

Kinetic  studies  were  conducted  on  both  single  and 
bisolute  systems  for  single  adsorbent  sizes  as  well  as  for  a 
mixture  of  adsorbent  sizes.  The  results  were  predicted  by  a 
computer  program  written  by  Mathews  (1975)  for  single  soiure 
system.  For  multicomponent  adsorption,  results  were  predicted 
by  a  program  written  by  Mathews  (1983).   It  was  seen  from  the 


predictions  that  Sauter  mean  diameter  gave  better  predictions 
than  geometric  mean  diameter  for  both  individual  size 
fractions  and  mixture  of  sizes. 

Fixed  bed  studies  were  conducted  for  single  solute  as 
well  as  for  bisolute  systems.  A  mathematical  model  of  fixed 
bed  was  modified  by  Mathews  (1984)  so  that  the  equilibrium 
could  be  represented  by  the  three  parameter  isotherm. 
The  model  was  employed  to  compare  the  performance  of  fixed 
beds  for  constant  temperature  and  pH  but  varying  particle 
sizes.  Breakthrough  curves  were  analyzed  for  geometric  mean 
diameters  and  for  Sauter  mean  diameters.  In  the  case  of 
mixture  of  sizes  the  model  was  analyzed  for  stratified  bed 
layers  also. 

It  was  seen  that  when  the  stratified  bed  layer  model  was 
used  the  column  predictions  were  better  than  for  either 
geometric  mean  diameters  or  5auter  mean  diameters.  However. 
it  was  seen  from  the  results  of  the  predictions  that  the 
representation  of  particle  size  distribution  by  the  Sauter 
mean  diameter  was  better  than  that  by  geometric  mean 
diameter . 


